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Abstract 

Strong long-scale gravitational waves can explain cosmic acceleration within the context 
of general relativity without resorting to the assumption of exotic forms of matter such 
as quintessence. The existence of these gravitational waves in sufficient strength to cause 
observed acceleration can be compatible with the cosmic microwave background under rea- 
sonable physical circumstances. An instance of the Bianchi IX cosmology is demonstrated 
which also explains the alignment of low-order multipoles observed in the CMB. The model 
requires a closed cosmology but is otherwise not strongly constrained. Recommendations 
are made for further observations to verify and better constrain the model. 
Keywords: Dark energy, acceleration, gravitational waves, Bianchi model, CMB, Axis of 
Evil, CMB cold spot. 
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Part I 



Introduction 



1 Background 



The observational confirmation that the universe has been expanding from a condition of extreme 
density and minute size since some point in the finite past represents a major triumph of Einstein's 
theory of gravitation in providing an elegant explanation for cosmology, without the addition 
of exotic, heretofore-unobserved substances or fundamental forces. This notion has however 
faced a serious challenge since Riess's 1998 discoverypQ of cosmic acceleration. The purpose of 
this research is to evaluate the question: can the back-reaction of long-wavelength gravitational 
waves in a closed universe contribute to cosmic acceleration while remaining compatible with 
observational constraints? 



1.1 Standard cosmology predicts an expanding universe 

The full Einstein equations reacQ 



where g^ is the metric tensor, R^ is the once-contracted Riemann tensor, R is the Ricci curva- 
ture scalar, T^ v is the energy-momentum tensor, A is the "cosmological constant" and the constant 
k = 8ttG/c 4 rj 2.08 x lCT^kg^m^s 2 . The Bianchi identity guarantees T xl = T 22 = T 33 so, in 
a Gaussian (g 00 = 1) and synchronous (g oi = 0) coordinate system we have(4]: 



(1) 



R°--R=kT°+A 



-R =kT>* - 2A 



(2) 
(3) 



lr Throughout this document, indices written with Greek letters fi, v etc. run over 0,1,2,3 and indices written 
in Roman letters i, j etc. run over 1,2,3. The sign of the metric tensor reads +,—,—,—. 
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Cosmological parameters 



In discussions of cosmology it is conventional to track the expansion of an isotropic metric by 
introducing a "scale factor", a positive function of time only. In general the scale factor has 
no specific geometric meaning other than to compare distances in the metric at different points 
in time. Furthermore, the scale factor loses unique meaning when the universe becomes non- 
isotropic. We will denote this scale factor function as a (t) in analogy with its definition in the 
Robertson- Walker metric, where it appears as [5] 



ds 2 = dt 2 - a 2 (t) 



(dx i y + (dx 2 y + (dx 3 y 



l + \K 



{dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 



(4) 



and the symbol K has the value in a flat universe, 1 in a closed universe, and -1 in an open 
universe. In this case the Einstein equations read3 



— (d 2 + K) =ke + A 



-6- =k(e + 3p) - 2A 
a 



(5) 
(6) 



where e denotes the energy density of matter described by the energy-momentum tensor and p 
denotes the pressure of matter described in that tensor. Let the Hubble parameter be defined by 
H = a I a and let Ho be the value of H at the present time (this is how we will generally use the 
subscript 0). Because the universe seems flat and dominated by ordinary matter over small scales, 
it is common to move terms arising from K to the right-hand side of the equation, where they act 
as elements of an "effective energy-momentum tensor", and to state contributors to cosmological 
expansion as dimensionless parameters Qi in comparison to the "critical density" £ cr itical' t na t 
is, the energy-density of ordinary matter required for the universe to be flat: ^e cr itical = so 

3H 2 = ke Q + A - 3| 
Q K + ft M + Mr + = 1 



where 



M 



n K = 



ke /3H 2 
A/3H 2 
-K/al 



(8) 
(9) 
(10) 



. multiple observations, most recently by WMAP, have confirmed that VIr <C Qm [HI and so to 
the limit of the precision with which these quantities can be evaluated Qk + &m + = lH 

2 A single dot denotes a derivative with respect to t; two dots denote a second derivative with respect to t. 

3 Chernin, in 1811 . elegantly derives a description of the scale factor in an open Friedmann cosmology which can 
be used when there is a significant amount of relativistic matter in a cold matter-dominated universe. Chernin's 
equation is easily generalized to the closed Friedmann universe. 
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1.2 Simple cosmology predicts a decelerating universe 

If the scale factor a measures a distance, it is reasonable to say by analogy that a can be compared 
to a velocity and a an acceleration. In isotropic cosmology we define the "deceleration parameter" 
Q InB 



y ~ a 2 dtH ( ' 



. Dividing ([6]) by ([7]) we easily obtain 

l fc(e + 3 P )-2A 1 
Q= 2* C + A-ir/ a ' = 2° M ^ A (12) 

. A flat universe with no cosmological constant must always decelerate. While the properties of 
so-called "dark matter" remain undetermined, the localisibility of dark matter's distribution and 
its slow motion implies it can be treated as w = dust. 

We can also immediately say that in a universe with no cosmological constant, acceleration is 
possible under the condition 

-l—(l-K/a 2 ke)<0 (13) 



1.3 Observations say the universe is accelerating 

Acceleration in and of itself is not a newcomer to cosmology. The de Sitter cosmology [HI], 
discovered in 1917, is driven solely by a cosmological constant and consequently has a constant 
deceleration parameter of Q = — 1. Bondi, Gold & Hoyle's "steady state" universe [83] similarly 
accelerates with Q = — 1, this value being associated with a universe whose expansion is driven 
solely by a field whose energy density is not dependent on the size of the universe. With the 
proposal of "big bang" nucleosynthesis |84j and the subsequent discovery of the cosmic microwave 
background [85 , consensus came to settle on the simplest matter-filled model, the Friedmann 
universe. [86J 

Throughout the 1990s, astronomical observations began to indicate that the matter energy den- 
sity of the universe was far below the critical density, leading some (for example [87] ) to propose 
the resurrection of the cosmological constant in order to preserve the observed near-flatness of 
space. 

In 1998, Riess et al. published an analysis[T| of the light from a small number of type la supernovae 
with 0.16 < z < 0.62 and concluded from this set that the recent universe is accelerating with 
Qo = — 1.0 ± 0.4. Further observations and analysis (see Part [II]) have also provided evidence 
that the universe has Q < 0. 

4 Q has been denned with a minus sign for historical reasons. Q > denotes a decelerating universe; Q < 
denotes an accelerating universe. We have avoided the more common notation q in favor of Q to avoid confusion 
when interpreting the source material. 
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While Riess et at did not exclude the possibility of a universe with K ^ 0, the assumption of a 
flat universe remains predominant throughout the field of cosmology as observations, both from 
supernova data and WMAP, have shown that the universe is, on observable scales, very close 
to flat - although it is impossible to distinguish between a universe that is genuinely flat, with 
Qk = and one with fix very close to but not equal to zero. 



2 Dark energy 



Since the discovery of acceleration, numerous explanations for the phenomenon have been pro- 
posed, all depending on an isotropic field creating additional, invisible energy. Turner and 
Huterer[6j introduce the term "dark energy", analogous to dark matter in the sense that dark 
energy does not interact electromagnetically with ordinary matter and has the property of an 
energy density, to term this additional energy, which appears to make up over 70% of the total 
energy content of the universe |16|. 

The assumption of a flat homogeneous cosmology demands that cosmic acceleration comes from 
a cosmological constant or a scalar field. Most scalar theories for explaining cosmic acceleration 
fall into two classes: an exotic form of matter with negative energy density, and surrender of the 
cosmological principle. Other scalar theories sacrifice different assumptions, such as homogeneity, 
or invoke more exotic explanations unsupported by laboratory physics. 



2.1 Cosmological constant 

The simplest most familiar variation on the Robertson- Walker cosmological model which allows 
an accelerating universe is the "ACDM" model - a universe dominated by "cold" (non-relativistic, 
p = 0) matter with both baryonic and dark components, and with the existence of a non-zero 
cosmological constant. In such a universe the Einstein equations read[5j 



SH 2 =ke + A (14) 

_6-=fe e -2A (15) 
a 

so when fee/A is small such that (he/A) 2 is negligible, that is, the universe is dominated by a 
cosmological constant, 

1 fee - 2A 3fce . s 

Q =2l7TA *" 1+ 2A < 16 > 

which at first glance appears to neatly explain Riesse et aUs result. However, as will be shown 
(see SECTION O, the case for a cosmological constant is not definite. Furthermore, the theoretical 
background explaining the strength of the cosmological constant is not well developed, relying 
on an understanding of quantum gravity which does not yet exist [72] • While the cosmological 
constant can always be said to have a "right to exist" in the Einstein equations, current physics 
does not explain why it should have any particular strength and as such the cosmological constant 
should be treated as the simplest form of a scalar field of exotic matter. 
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2.2 Quintessence 

More general than the cosmological constant but similar in structure is the proposal of "quintessence" [BJ, 
a novel form of matter with a time-dependent equation of state that can take on negative val- 
ues. Many forms of these have been proposed; one form of these, for example, is the "Chaplygin 
gas" [55], which has equation of state p = — A/e for A > 0. Quintessence theories are particularly 
motivated by the idea that acceleration is a cosmologically recent phenomenon, noting limited 
data (see Part[ITJ that the equation of state of dark energy may be evolving with time. 

Any formulation of quintessence must be regarded as highly speculative. At the most fundamental 
level, all theories of quintessence propose the existence of a kind of matter which: 

• has never been observed experimentally; 

• does not interact with ordinary matter via the electromagnetic force; 

• has a negative equation of state, that is, a positive energy density produces a negative 
pressure; 

• plays a prominent role at current energy levels, as opposed to exotic effects (e.g. unification 
of forces) thought to have taken place only in the very early universe. 

In the absence of any compelling experimental evidence whatsoever for any kind of quintessence, 
quintessence and quintessence-like models should be rejected as definitive explanations for dark 
energy. 



2.3 Local inhomogeneity 

A more mundane explanation which has been offered for acceleration is the "Hubble bubble" [Tll61|. 
regions of lower density in the intergalactic medium. If the vicinity of the Milky Way had lower 
matter energy density, expansion in its vicinity would increase[93j, causing the illusion of cosmic 
acceleration. 

Not only would the density deficit in such a "bubble" have to be quite large in order to cause 
acceleration, but the theory, which has the advantage of requiring no new physics, supposes 
either the existence of a rare or unique void that the Milky Way happens to be in - a violation 
of the cosmological principle in the sense that it makes observers in the Milky Way privileged - 
or a preponderance of voids whose presence makes the universe inhomogeneous not just in small 
patches but on average. [511 IM] 



2.4 Exotic models 
Modified relativity 

Some proposals to explain dark energy propose modifications to the Einstein equations. The best- 
known of these is the Cardassian Expansion model[95j, which proposes time-dependent variation 
of the equation of state of matter. The Cardassian model is of particular interest in that it 
proposes an equation for the density perturbation 
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S 2 K = 



(17) 



for unknown constant s, which equation begins to resemble that for weak gravitational waves 
in a closed universe (cf. equation (|153|l ). Like Chaplygin gas and the "DGP" model[§S], the 
Cardassian model justifies itself based on theories about higher-dimensional manifolds which 
remain untested. 



Topological defects 

The existence of cosmic strings would change the overall equation of state of the matter in the 
universe by a constant [571 \§7\, creating acceleration through simple deviation from the Friedmann 
model. While theories of cosmic inflation predict the formation of cosmic strings and other 
topological defects, such defects remain completely undetected. 



3 Tensorial theories for acceleration in a flat universe 

If we wish to preserve the theory of relativity and cosmic homogeneity, while at the same time 
relying only on effects with good experimental basis, scalar fields appear to be excluded as an 
explanation for acceleration. Ergo within the context of general relativity the next place to search 
for an answer is in tensor theories, which include the possibility of gravitational waves. 

Lifshitz's theory of cosmological perturbations [79J appears to exclude tensorial answers to the 
problem of acceleration: gravitational waves have the same equation of state as radiation, and 
local clumps gravitational waves in the theory (where "local" means bounded within an region 
smaller than the radius of curvature of the universe) both decay rapidly and collapse spatially. 
Rodrigues[113] takes a first step in discussing anisotropic dark energy, but limits his analysis to 
a flat universe and thus creates the problem of an anisotropic "big rip". 

A high-frequency gravitational wave background has been proposed[92| as the source of cosmic 
acceleration. While the authors' analysis appears initially promising, similar to many scalar dark 
energy candidates the theory relies on the existence of an inflation-induced gravitational wave 
background that remains only hypothetical. Furthermore, the authors obtain their result by 
selection of an averaging scheme without mathematical rigor - surely choosing a mathematical 
model based on the desired results cannot be considered physics. At any rate, the strength of 
the background inflationary theory predicts is not sufficiently great to explain the observed large 
acceleration. 
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4 Introduction 



The theory of tests to evaluate the deceleration parameter using supernovae as standard can- 
dles began with Wagoner |45| in 1977. Starting from assumptions of an isotropic Friedmannian 
cosmology which is not necessarily flat, Wagoner notes the approximate relation 



±(l + Q )z 2 + O(z 3 ) 



(18) 



which, when Ho and z are known, relates the deceleration parameter to the distance <1e as 
determined by the dimming of the supernova (where Wagoner was originally discussing Type II 
supernova events jf[ This relation is valid when z is small such that z 3 is negligible, limiting its 
usefulness above z ~ 1, and requires the assumption of only small changes in the Hubble constant 
Hq (that is, in a Friedmann cosmology apjaF evaluated near the observer) on the interval from 
z « to z w 1. 

Type la supernovae are thought to be a "standard candle" for the measurement of distance 
and redshift; that is, supernovae of that type are thought to possess spectral and luminosity 
curves which are nearly identical. Therefore, observation of extragalactic type la supernovae is 
believed to produce reliable information on both the distance of the event (noting that brightness 
diminishes as the inverse square of distance) and the redshift of the distance associated with the 
event (through the change in the peak of the supernovae's spectra), with redshift z related to 
the scale factor ap by 

z _|_ ^ _ ap (^observation) Qg\ 
aF (^emission) 

. Analysis of a statistically unbiased dataset of z (£) therefore gives empirical information on 
Q(t). ' 

Colgate [U] proposed that Type I supernovae should be used to measure the deceleration param- 
eter in preference to Type II supernovae. Type I supernovae, specifically the "Type la" whose 
mechanism is thought to be the accretion of matter onto the surface of a white dwarf star, are 
understood to have a well-defined typical absolute magnitude and spectrum, and assuming this 
is true the distance to and redshift of a given Type la supernova event (SNe) can easily be de- 
termined by fitting its light curve to standard templates. Therefore, with a sufficient sample of 
extragalactic supernovae of z < 1, the parameters H and Q can be measured directly. When an 
isotropic cosmology with constant deceleration parameter Q — Qo is assumed, knowledge of Hq 
and Qo are sufficient to typify the parameters of the universe |4J. 

With the advent of modern optical astronomy such as adaptive optics [46J and space-based optical 
telescopy [47], such surveys have become possible, but have produced results contradicting the 
standard, cold matter-filled Friedmann model of cosmology. 



^Equation fl 18D is of course a generalization of the famous distance-redshift approximation HqcLe ~ z.[H 
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5 Surveys of acceleration 

Cosmological studies measuring Q have been ongoing since 1997 and consist of analysis of 
redshiftsHJ 121 IS1 [III IISl IISl [T^l BSl ISSl ISSI of type Ia supernovae. 

The High- z Supernova Search Team's initial study of the deceleration parameter[T] was the 
first large study to call attention to the problem of acceleration. Working from a sample of 
sixteen supernovae (four of which were well-observed "high confidence" sources), the most distant 
with z = 0.97, Riess concluded that the universe has Qq < to high confidence, although 
the measurement of Qq itself possessed a high degree of uncertainty. Riess also noted the high 
sensitivity of the result to individual data points. Oddly, the authors dismiss the closed cosmology 
despite the data indicating it as preferred[TJ Fig. 7]; however the size of their experimental error 
precludes real evaluation of spatial curvature. 

The Supernova Cosmology Project (SCP) made an earlier attempt to evaluate Qo with the use 
of supernovae [55]. This small survey (n — 7) on relatively nearby supernovae found a result 
inconsistent with those that followed it, giving results consistent with a universe with no dark 
energy and with too high a degree of error to meaningfully evaluate the geometry of the universe. 

In contrast, the Supernova Cosmology Project's 1998 evaluation (151 EE9"] of 42 Type Ia SNe added 
further evidence that the universe was accelerating, and also makes note of the surprising coin- 
cidence of the energy density fl\ 's near-equivalence with the total energy density in the current 
epoch. The SCP also failed to consider the closed cosmology despite supernova data favoring 
it 021 Fig. 7]. 

The ESSENCE[TT] survey was expressly designed to examine cosmic acceleration and detected 
102 type-la supernovae from 0.10 < z < 0.78, of which 60 were used for cosmological analysis. The 
initial analysis of ESSENCE assumed flatness of the universe. ESSENCE'S observational fields 
were deliberately chosen to overlap the areas of previous surveys and to lie within ten degrees 
of the celestial equator; all were also between 23:25 and 02:33 Right Ascension. Combining data 
from ESSENCE, SNLS and other sources[52j led to a conclusion consistent with other analyses. 
Exploration of more exotic models [53J found that no model of those tested was a good fit for 
ESSENCE'S data. 

The Supernova Legacy Survey (SNLS)[12J recorded 472 type-la supernovae. While analysis of the 
SNLS dataset[T3] provides results consistent with a universe driven by cosmological constant, the 
uncertainty on analysis of a time-dependent component to the equation of state of dark energy 
is very large; their analysis also does not consider a closed universe as a possible model[50]. 
Furthermore, the SNLS team also note the presence of two outliers and only 125 of 472 events 
were used to evaluate cosmology. SNLS observed SNe in four fields, one of which (field 3) is far 
above the plane of the celestial equator at 52 degrees declination; this and [53] 's northern field are 
the only fields with multiple observations in a small area more than 20 degrees from the celestial 
equator surveyed to date. SNLS also notes[50l section 5.4] that the values of evaluated in 
the four fields are compatible only at a 37% confidence level - a surprising result given that each 
SNLS field contains at least 60 SNe in quite small (one square degree) areas. 

The Hubble Space Telescope or HST survey of supernovae, published in 2004 and reviewed by 
the Supernova Cosmology Project[T3] observed twenty type-la supernovae with redshifts 0.63 < 
z < 1.42. While the number of SNe observed is small, the HST survey has the advantage of 
covering a wider area of sky than other SNe surveys. Analysis of the HST dataset suggests a 
rapidly-evolving dark energy field, although with very high error on measurements greater than 
z = 1 due to the small (n — 10) sample size it is impossible to take these results as anything 
more than suggestive. HST slightly favored a closed model of the universe, when considering 
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interpretations of data that allowed Qk =h 0- 

The Supernova Cosmology Project's 2008 analysis of supernova data[51J made a analysis of 
combined SNLS, ESSENCE and HST data, and attempted to analyze the data in the context 
of a theory of a time-dependent equation of state for dark energy but concluded "present SN 
data sets do not have the sensitivity to answer the questions of whether dark energy persists to 
z > 1, or whether it had negative pressure then." The analysis rejected 10% of all SNe from the 
combined data sets as outliers, many based on their failure to fit with a nearby Hq; Kowalski 
et al. 's rejection of outliers also shifts their analysis from one favoring a closed universe to one 
favoring a flat one [51, Fig. 11]. 

Further work by Riess et. al. |54l [55] produced the so-called "gold" dataset of SNe, a group of 
supernova events with particularly clear light curves with 33 at z > 1. These supernovae were 
observed in two small (one square degree) fields. |54] claims a great reduction in the uncertainty 
of the Hubble parameter at z > 1 but the Hubble parameter measured in the extended "gold" 
set gives a value for the Hubble parameter not reconcilable with that in the [56] dataset. Riess 
et. al. conclude that w is negative (with large experimental error) in the region 1 < z < 2, then 
attempt to extrapolate the behavior of dark energy back to z — 1089. 

Sollerman et al 's analysis of the Sloan Digital Sky Survey-II supernova datapU [5] is the most 
recent analysis indicating cosmic acceleration. SDSS-II observed 103 type-la supernovae in a 
long, narrow strip along the celestial equator, including many from lower redshifts than had been 
previously examined in detail; Sollerman et. al. also made use of data from the HST, ESSENCE 
and SNLS surveys, bringing the total number of SNe examined to 288. The primary conclusion 
to be drawn from SDSS is the sensitivity of cosmological measurements to the specific analysis 
technique used[60 ; analysis of the data with two different curve-fitting algorithms produce two 
different, albeit somewhat compatible, results. 

Further obscuring the neatness of measuring Q, Jha et. al. notcd[59[ that the uneven local 
distribution of galaxies, specifically the existence of voids, can lead to a mis-estimation of Hq on 
the order of 6.5% for a given galaxy. 

Finally, of note is the WiggleZ dark energy survey (15] 143] . WiggleZ is the most extensive redshift 
survey thus far conducted, with some 280,000 galaxies with 0.2 < z < 1.0 used as sources. 
WiggleZ also covers a wider area of sky than previous surveys, examining some 1000 square 
degrees in multiple windows around the sky. Two of WiggleZ's windows overlap with SDSS-IFs 
survey area, so while WiggleZ is ongoing preliminary results [57, 58] can be used to improve the 
evaluation of Q by improving precision on measurements of z of SNe host galaxies. The authors 
of [58] note that "the redshift-space clustering pattern is not isotropic in the true cosmological 
model", attributing the variation to "the coherent, bulk flows of galaxies toward clusters and 
superclusters". Analysis by the WiggleZ team of pre-existing SNe datasets, using the new, more 
precise data on galaxy redshifts they obtained, reconfirms the fact of acceleration, and generates 
results consistent with other surveys, but the data lack sufficient precision to determine the 
history of Q. 

Table [3] details the sky locations of SNe and galaxies used in the determination of acceleration; 
Figure [3] presents these locations graphically. Table [T] in Appendix summarizes the results of 
these surveys. 
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Survey 



Supernova Cosmology Project 1997[55j 



High-z Supernova Search Team[l 



Supernova Cosmology Project 1998[49J 



HSTI14I 



ESSENCE! 11 



Supernova Legacy Survey [121 150| 



ESSENCE + SNLS[52) 



Supernova Cosmology Project combined[51 



Riess "gold" sample [561 154] 



WiggleZ [58 



No of SNc 



7 



16 



42 



20 



102 



125 



162 



307 



16 



557 



0.35 < z < 0.46 



0.16 < z < 0.97 



0.18 < z < 0.86 



0.63 < z < 1.42 
0.10 < z < 0.78 



0.015 < z < 1 



0.015 < z < 1 



0.015 < z < 1 



1.25 < z < 2 



O.K z < 0.9 



Survey 


tl\ 




tlx 


Supernova Cosmology Project 1997 


0.06t££ 


u - y *-0.28 


(dne) 


High-z Supernova Search Team 


72 +u - T2 

u - ' z -0.48 


u - z ^-0.24 


(dne) 


Supernova Cosmology Project 1998 


n 79 +0.08 
u - ,z -0.09 


o.28±8:g 


(dne) 


HST 


n 71 C+U.U36 
u -' lo -0.057 


0.286irol 


-0 001+ tu)3V 


ESSENCE 


726 +,UM) 


0.2741S 


(dne) 


Supernova Legacy Survey 


0.751 ±0.080 


0.271 ±0.020 


(dne) 


ESSENCE + SNLS 


n 7 oo+0. 018 
u - '°°-0.028 


267 +U.U28 
u - z0 ' -0.018 


(dne) 


Supernova Cosmology Project combined 


u. <eso_ 045 


285 +U U3U 

U.ZOO_Q Q3Q 


-0 010 +u uib 


Riess "gold" sample 


7l+ l) - l)3 


29 +u ub 

u - zy -0.03 


(dne) 


WiggleZ 


0.71 ± 0.03 


0.29 ±0.03 


(dne) 


Sloane Digital Sky Survey-II 


0.693 ± .042 1 
0.744 ± .041 2 
0.74 ± .02 3 


0.307 ± .042 1 
0.256 ± .041 2 
0.25 ± .02 3 


.04 ± .04 



Survey 


n hatt 
Vo 


WXO 


WXa 


Supernova Cosmology Project 1997 


u - 41 -0.17 


(dne) 


(dne) 


High-z Supernova Search Team 


-1.0 ±0.4 


(dne) 


(dne) 


Supernova Cosmology Project 1998 


-0.58+^ 


(dne) 


(dne) 


HST 


-0.5721^ 


—0 QQ7+ u;2bb 
u.jy 1 _ 293 


0.13±if 7 


ESSENCE 


-0 589 +u ' u ^ 


J-U*< _0.124 


(dne) 


Supernova Legacy Survey 


-0.620 ±0.060 


-1.023 ±0.087 


(dne) 


ESSENCE + SNLS 


-o.eoolffi 




(dne) 


Supernova Cosmology Project combined 


-o.642±8:SSf 


-i.ooi^;^ 


(dne) 


Riess "gold" sample 


u - JD -0.02 


-1 02 +lU3 




WiggleZ 


-0.56 ±0.02 


(dne) 


(dne) 


Sloane Digital Sky Survey-II 


-0.539 ± .021 1 
-0.616 ± .021 2 
-0.61 ± .01 3 


-0.76 ± .18 1 
-0.96 ± .18 2 


(dne) 



Table 1: Summary of results from surveys indicating acceleration 
"dne" = "Does not evaluate". *: J56] / attempts to analyze w a with several different constraints 
but provides no numerical figure for its estimate of w a 's value, f : Where not explicitly stated in 

the source, t}\ is evaluated from Qm + tlx + = L t- <3o^ = htljtf — tl\. (1): MLCS2K2 
evaluation. (2): SALT -11 evaluation. (3): ACDM model evaluation. 
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Figure 1: Sky positions of supernovae used as evidence for acceleration 
Surveys of cosmic acceleration cover a limited portion of the sky, and data are divided into two 
contiguous, antipodal regions. Most data has been collected in a small area of the sky near the 
equator. Triangles: Riess 1998 supernovae. Five-pointed stars: HST SNe. Six-pointed stars: 
Riess "gold" dataset. The long, thin strip centered on 0,0 is the SDSS-II survey area. Other 
boxes are the SNLS and ESSENCE survey areas. 

6 Analysis 

Analysis of supernova data is, in one sense, quite consistent: all surveys apart from |55| agree that 
for z < 1 we have a deceleration parameter Qo — —0.6. Deeper analysis suffers from a lack of data 
at high redshifts and large numbers of free parameters in cosmological models, especially when 
more exotic models are considered. Meanwhile, while most surveys indicate that the acceleration 
in recent times acts as though driven by a cosmological constant, with an equation of state 
compatible with wx = — 1, the results from [60 show that this can be the result of the prior 
assumptions made about the model of dark energy. 

No definitive statement can be made about the evolution of Qq over time from the information 
thus far available, particularly not statements connecting the state of cosmic acceleration now 
with the state of acceleration at the epoch of last scattering. 

Nor can any definitive statement be made about cosmological models, other than to say that the 
most conservative, ACDM model fits the data at best inconsistently. Few studies of supernova 
data on acceleration examine the question of curvature in depth. 

The majority of SNe data is collected from a single patch of sky: the field bounded by RA 22:00, 
RA 04:00, Dec +1°15' and Dec -10°00' (the "highly-observed field"). This area comprises 1350 
square degrees, or only 2.1% of the sky. Surveys taken in small fields outside the highly-observed 
field, such as the Riess "gold" dataset, have high internal consistency, while surveys covering larger 
areas of sky have much lower consistency; the "gold" dataset contains the same number of SNe 
as the pQ but has a standard error less than a tenth the size. It is also telling that the four SNLS 
fields produced results that correlated poorly (37% confidence) with one another [S"U]. where two of 
the SNLS survey regions are well outside the highly-observed field. Compounding cosmographic 
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bias, many of the remaining SNe observations are located in a region of sky antipodal from the 
highly-observed field; any vector or tensor contribution to cosmic dynamics will be dominated 
by dipole and quadrupole terms, and as such be seen with equal or opposite magnitude in the 
antipodal direction (that is, if we observe a change in Q of AQ along the redirection, we should 
expect a change of — AQ in the event of a vector contribution, or AQ in the event of a tensor 
contribution, along the — x l direction). 

There is, furthermore, no SNe data whatsoever from above Dec +62°or below Dec —37°. The 
authors of [58] note a variation in the apparent Hubble parameter for galaxies in this equatorial 
band (no WiggleZ region lies further north than Dec +8° or Dec —19°); variation to the Hubble 
flow could potentially be even greater outside this region. There is also no evaluation of whether 
the Hubble flow remains isotropic beyond z = 0.3 |77] . 

Indeed, Zehavi et. at comment[5H] on the lack of sky coverage in their analysis of local Hubble 
flows, noting that "sparse sampling and the incomplete sky coverage (especially at low Galactic 
latitudes) may introduce a bias in the peculiar monopole due to its covariance with higher multi- 
poles". While the fact of greater redshift in the range where acceleration can be measured should 
overcome the peculiar velocities of galaxies, the data problem remains. 



7 Conclusions 

Many reasonable constraints prevent a full-sky survey of supernovae. In the optical band, much 
of the sky is obscured by the "zone of avoidance" created by the plane of our own galaxy [53]. 
The so-called "Great Attractor", certain to be a region of particularly high peculiar velocities and 
therefore great shifts in the apparent Hubble parameter, lies in this zone [64]. Furthermore, with 
only a single space-based optical observatory (the Hubble Space Telescope) operating, detailed 
observation of the sky is restricted to those latitudes accessible by ground-based observatories, 
none of which are located in Arctic latitudes. However, the directional deficit of SNe surveys, 
aggregated together, cannot be ignored. 

In the light of Tegmark et. al.'s discovery [28] of a preferred axis to the CMB quadrupole, and 
Land & Maguiejo's subsequent observation|32] of a preferred axis in higher multipole moments 
aligned with the the quadrupole (the so-called "Axis of Evil") , the default assumption should be 
that anisotropic acceleration is not ruled out. Indeed, the prominent CMB "Cold Spot" [34] falls 
within the highly-observed field, although no surveys or SNe are located exactly in its direction. 

As such, Wagoner's assertion of the cosmological principle as "statistically valid" [45] has been 
misapplied by analysts of acceleration data. A tensorial theory of cosmic acceleration would pre- 
serve homogeneity, in the sense that every observer sees "the same version of cosmic history" [2"T] . 
at the expense of isotropy in the form of spherical symmetry. 

More fundamentally, all studies of cosmic acceleration to date operate on the assumption that 
acceleration is isotropic, that is, that the acceleration field is equal in every direction, and therefore 
must be explained either by a cosmological constant or a scalar field. As Mortsell and Clarkson 
note, "[a]t best this gives a small error to all our considerations; at worst, many of our conclusions 
might be wrong" [61]. In particular, the data as presented cannot distinguish between a scalar- 
field theory of acceleration, a vector-field theory of acceleration, a cosmological constant theory 
of acceleration, and a time-dependent tensor-field theory of acceleration. 

Meanwhile, the simplest theory of acceleration, a cosmological constant, is challenged on two 
fronts: not only is fl^s value far out of line with that predicted by theory]72], but while its 
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equation of state is close to wx = — 1 measurements have tended to favor a value slightly smaller 
than -1. 

It is interesting to note that when £Ik is evaluated, supernova data favor a closed universe 
(although always in a manner compatible with a flat universe) ; this conclusion is consistent with 
the curvature parameter evaluated by WMAP[62J. 



7.1 Recommendations 

In light of these weaknesses of the current information on cosmic acceleration, the following 
program is recommended: 

The data already in hand should be re-evaluated to look for signs of angular dependence in the 
Hubble parameter. In particular, the SNLS fields 2 and 3 and the "gold" sample fields should 
be examined in and of themselves in order to build up a map of H as a function of both z and 
direction in the sky. The rejection of certain SNe in [51] should be re-evaluated in light of possible 
inadvertent obscuring of evidence for angular dependence in H. 

Analyses of SNe data should always consider the possibility of a closed or open universe as well 
as a fiat one. 

Additional SNe surveys for redshifts .15 < z < 2 should be carried out in unexamined areas of sky 
not obscured by the plane of the galaxy, such as for example the celestial north and south poles. 
The optimal region for these surveys is in rings located 90° from the center of the highly-observed 
field, which will maximize the difference in the event of a tensor-field (that is, gravitational-wave) 
acceleration. 

In light of this need and the lack of ground-based observatories, as well as the Zone of Avoidance, 
priority should be given to the Wide Field Infrared Survey Telescope ("WFIRST") project [65]. 
which incorporates the Super Nova/ Acceleration Probe[66, 67J and Joint Dark Energy Mission[68, 
161?] . This telescope is currently scheduled to be launched in 2016. 

As WiggleZ continues, its data on galactic redshifts should be examined for angular dependence 
as well. The completion of WiggleZ will provide invaluable information on baryon acoustic 
oscillations which will make possible the charting of the history of H and Q at much higher 
redshifts than is possible through the examination of supernova data. 

Zhao et al. have also noted the possibility of using the Einstein telescope as an instrument for 
examining dark energy through the use of gravitational wave emissions from colliding binary 
objects as a "standard siren" analogous to the standard candle of type la SNe |lllj . The largest 
binary systems, in the period prior to merger, could be parameterized by radio telescopy [116J. 

Cooray and Caldwell[102j, implicitly identifying the same problem of lack of angular coverage as 
we note herein, propose a program of near-redshift surveys covering a large but practical area of 
sky which could also provide the relevant information with existing facilities. 

Overall, the need is underscored for new theories of acceleration, particularly ones that attempt 
to explain acceleration through the action of tensor perturbations in a closed universe. Wagoner's 
formula (equation [18]) and its generalizations must be generalized further, to take into account 
the possibility of anisotropic fields as the cause of anisotropic cosmic acceleration. 
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Part III 

The Bianchi IX cosmology 

In pursuit of a theory within the context of unmodified general relativity which can explain 
cosmic acceleration while remaining compatible with the cosmic microwave background, we wish 
to relax as few constraints on our cosmological model as necessary. Therefore while having 
sacrificed the requirement of isotropy in the sense of spherical symmetry in the dark energy field, 
we wish to retain a stronger^ ss. 116] condition of the Copernican principle on our space, that 
of homogeneity |H Chap. 13 sec. 1]. It is also desirable to have a model whose limiting case is 
a Friedmann cosmology, in order to explain the almost-isotropic (that is, almost-Friedmannian) 
character of the CMB. Furthermore, a model which is spatially closed, in order to match models 
favored by CMB and SNe data, is desirable; such a model would, if complying with all other 
conditions, have a flat universe as a limiting case in the limit of an infinitely large radius of 
curvature. 

Bianchi showed [70] that there exists exactly one homogeneou^l space with a closed Friedmann 
universe as a limiting case: the Bianchi type IX cosmology, for short "Bianchi IX". 



8 The Bianchi classification scheme 

Bianchi observed that all three-dimensional homogeneous spaces could be classified into nine 
types, based on categorization of the symmetries, that is the Killing field, in each space. Behr 
noted[71j that this categorization scheme could be simplified to filling a parameter space of four 
parameters: one running over the real numbers and three reducible to the sign function sgn(x). 

Consider some space with metric ds 2 — dt 2 — gijdx 1 dx 3 (that is, a space in Gaussian coordinates) 
where g.^ = gij (t, a; 1 ) . If the sub-space with metric tensor g^j is homogeneous, then there 
exists a set of vectors that solve ^.j + = 0; these are the Killing vectors of the space[4]. 
In a homogeneous space, these vectors will (where [ , ] is a commutator) obey the commutation 
relationship 

[£i,Z j ] = titi-(jti = C k ij tk (20) 

where in an homogeneous space, the object is a constant pseudo-tensor, the "structure con- 
stants" of an homogeneous space, with the antisymmetry property C\n = C k j ss. 116]. 

We always have the freedom to perform separation of variables the functions gij ; let us do so by 
defining the matrix "f a b such that 

g ij (t,x k )=- lab {t)e ( ? ) (x k )ef ) {x k ) (21) 



6 A homogeneous space is a space such that for any two points in the space, there exists a geodesic, not 
necessarily of finite length, connecting those two points. 
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The 3x3 matrix e,- a) (x k ) is a triad[2l I112( ss. 980 of vectors ("frame vectors") which solve 

& (e\ a) ef Wda/) = (22) 

(in the language of linear algebra, the quantity e^dx 1 is a one- form on a homogeneous space). 

Furthermore, define the matrix such that e5 Q N = = <5*; from this it follows that e^ej = 
. From these relationships we can transform between any tensor and its decomposition into 
triads by saying that for some tensor A 1 ^ 2 ^ " l p , 



a)l{a.) 2 {a) 3 ...(a) m ( h l2 J 3 i \ ( {b) 1 (6) 2 (6) 3 J b )A /oo\ 

6) 1 (6) 2 (6) 3 ...(6)„ V e (a) i e (a) 2 e (a) 3 ■•■ e (a)„J ^! 6 'n e j3 ' ' " ^ ^ 



; therefore in an homogeneous space we can perform separation of variables on the partial differ- 
ential equations of general relativity and solve the time-dependent parts as ordinary differential 
equations. 

The frame vectors obey the properties 



eg " eg = C% c efef (24) 

[TJ. The structure constants C bc typify a homogeneous space and are given by the following 
rule fTTfl : 

C a bc = e bcd n ad + 6 d c a b - S d a c (25) 

where the object n ah is a diagonal matrix diag (n^ 1 ', n^ 2 \ n( 3 )) and a a is the vector (a, 0, 0), the 
values of this matrix and vector typifyied by the underlying cosmology (Table [5]) . 

The cosmologies of Bianchi types I, V, VIIo, VII a and IX are of particular interest as they have 
isotropic spaces as limiting cases; specifically, a universe with metric 

ds 2 = dt 2 - a^e^ef (26) 

is a flat K = universe for Bianchi type I, an open K = — 1 universe for Bianchi types V, 
VIIo, VII a , and a closed K = 1 universe for Bianchi type IX[2J QI3 HI]- Bianchi IX is the only 
homogeneous closed cosmological model in the context of general relativity[104j. 



7 Indices from the beginning of the Latin alphabet (a, b, c,...) denote triad indices; indices from the middle of the 
alphabet (i, j, k,...) denote regular indices. Where the two are mixed or the application is otherwise ambiguous, 
triad indices are enclosed in parentheses; in this work, this notation never means the tensor symmetrization 
operation. 

8 The widespread Fourth Revised English Edition of [2] contains numerous serious typographical errors in the 
section introducing the tetrad formalism. The Russian-language Seventh Corrected Edition \ll2\ contains the 
correct formulas. 

9 The symbol e a f, c represents the Levi-Civita symbol defined such that £123 = 1 
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Bianchi type 


a 


n y ' 


n K 1 


^3 
TV- ' 


I 














II 





1 








III 


1 





1 


-1 


IV 


1 








1 


V 


1 











VIo 








1 


-1 


VI a 


a 





1 


-1 


VII 





1 


1 





VII a 


a 


1 


1 





VIII 





1 


1 


-1 


IX 





1 


1 


1 



Table 2: The Bianchi classification scheme 
Constants for the different homogeneous spaces of the Bianchi classification scheme [21 |2"T1 ITT1 1104] . 
The quantity a runs over the real numbers. This parametrization is not unique (we could, for 
example, have chosen (—1, —1, —1) for (n^ 1 ', rv- 2 \ nA 3 )) in the type IX space). 



9 The Kasner universe 



In order to illustrate the possible effects of an anisotropic but homogeneous cosmology on cosmic 
dynamics, we will consider a Bianchi type I cosmology that generalizes the Friedmann cosmology: 
the Kasner uni verse |76|: [H ss. 117]. 

Let our metric read 



ds 2 = dt 2 - 1 2 ^ (dx 1 ) 2 - 1 2 ^ (dx 2 ) 2 - 1 2 ^ (dx 3 y 



(27) 



where pi , p2 , P3 are constants. In a co-moving coordinate system we quickly arrive at the following 
set of Einstein equations: 



[(pi +P2+P3) - 


(pi- 


VP\ 


+Pi)]t~ 2 


=-k(e 
2 V 


f 3p) 


(Pl- 


VP2 - 


Vp-i 


-I) Pit- 2 


— — k (p 
2 yy 


-e) 


(pi - 


VP2 - 


Vp-i 


-l)P2^ 2 


= — k (p 
2 yP 


-e) 


(pi - 


VP2 - 


Vp-i 


-I) Pit- 2 


= — k (p 
2 yP 


-e) 



(28) 
(29) 
(30) 
(31) 



. These equations necessitate either an isotropic but unphysical (p — e) universe or a vacuum 
(e = p = 0) universe, in which we have either the trivial solution p\ = p 2 = P3 = (Minkowski 
space) or the more interesting solution 



P1+P2+ P3 =P 2 i+pI+pI = 1 



(32) 



This solution admits a parametrization of p\ 1 p 2 ,P3 such that (if we choose p\ < p 2 < P3) 
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Pi = — u/ (l + u + u 2 ) 

P2 =(1 + U) I (l+U + U 2 ) 
p 3 =U (1 + U) I (l + U + U 2 ) 



(33) 



where u > 0; these relations have the nice symmetry property that pi (u) = pi(\/u). 



9.1 Scale factor 



The scale factor a does not necessarily have an intrinsic meaning, but instead compares distances 
as a function of time. In an isotropic cosmology such as the Friedmann model a can be given a 
real geometric meaning; in an open or closed Friedmann universe, the scale factor appears simply 
in the Ricci curvature of space R l - = (2if/a 2 ) <5* and as such can be regarded as the radius of 
curvature of the universe. In particular, in a closed isotropic universe a can be considered to 
have the direct physical meaning of the radius of curvature of the spherical space, so in a closed 
universe one could meaningfully say "the radius of the universe is a". 

When space is no longer isotropic, the definition of scale factor breaks down. It is, of course, 
possible to define any positive function as "the" scale factor. Grishchuk et a/.|10|, section 4], for 
example, use a metric 



7n 

722 
733 



-a 2 e 2a 



a 2 e 2 ? 



2 2t 

a e ' 



(34) 



and propose the definition 



12 



labV 



ab 



(35) 



in the context of a vacuum cosmology, motivated by the coincidence of this definition of the scale 
factor with one the authors introduce in separating the Bianchi IX metric into background and 
gravitational-wave parts. The authors also discuss a definition of scale factor such that 



(det7 afc ) 1/3 



(36) 



. This definition has the advantage that it relates the scale factor to a definite physical quantity, 
a volume element, but it contains a deeper flaw: with such a definition in place the Einstein 
equations admit no solution other than the background solution at quadratic and higher orders. 
If we define the quantity 



5 = a + /3 + ~/ 



(37) 
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then 



(det7a 6 ) 1/3 



(38) 



. In either case, though, discussion of definitions of a attempt to solve a problem that does not 
exist. The question of what definition of scale factor to select is analogous to the question of 
which of the orthocenter, incenter or circumcenter of a triangle is the "true" center. Consequently, 
attempting to extract a single scale factor - and thus a single Hubble parameter or a single 
deceleration parameter - from anisotropic Einstein equations is a fool's errand. 

We can, if we wish, split the metric (fM]) into isotropic and anisotropic parts by noting that the 
quantity aj?e* is isotropic and that any two of the quantities a — (3, a — 7 and (3 — 7 combined with 
a^e* contain all the information needed to describe the metric [22]; pursuing this route would be 
a distraction from our main task. 

Instead, let the notion of scale factor a, Hubble parameter H and deceleration parameter Q be 
generalized. In a homogeneous cosmology with a diagonal metric, define the following matrices: 
the generalized scale factor, 



( (7n) 1/2 


V 





(722) 







(733) 



(39) 



1/2 



(recalling that non-integer powers of a matrix are not defined, so we could not simply say 



a a b 



(lab) 



1/2n 



In a Bianchi I cosmology only, from this definition we can then define the redshift 



matrix (in homogeneous cosmologies other than Bianchi I the geodesic equations are non-linear; 
see Part [TV]) : 



z b a EEa ac ( m )a bc ( VE )~S b a 









n22(*fl) 
a22(tB) 








a33(*,R.) 



(40) 



where the subscript R denotes the function evaluated at the time of observation of light, and 
E denotes the function evaluated at the time of emission; and finally the generalized Hubble 
parameter and deceleration parameter: 



1 d 



an/an 
022/022 
033/033 



(41) 



Q» a EE jB^bc -K = - 



a n a n / (dn)" 







022022/ (022)^ 







033033/ (033) 2 



(42) 



. This approach is essentially a generalization of that developed by Barrow in [22]; the object 
(|4ip is closely related to the shear tensor [2"Tl 126] . The practical purpose of these definitions is to 
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provide a mathematical description of observed quantities; let e l be a unit vector pointing in the 
direction of observation. Then the redshift observed in the e l direction is given by 



z{e\t)^z a b e\ b) e { J a) ^e l (43) 

and similarly for other functions of the scale factor. Each of these functions can be averaged over 
the whole sky, these averages denoted by a bar: 

a . Jv J e3ds = -a ab ^ - 3 (-n + a 22 + a 33 ) (44) 

etc. 



9.2 Dynamics in the Kasner universe 



An observer in a Kasner universe will see the consequences of that universe's evolution. Exami- 
nation of the observational consequences of the Kasner universe provides an illustrative example 
of potential consequences of anisotropy in other cosmologies. 



Expansion 



Misner, Thorne & Wheeler argue[76j that the Kasner universe is expanding, as the volume element 
is always increasing: 



^ = j t ^~\\dx 1 dx 2 dx i = j t (t Pl+P2+P3 ) dx 1 dx 2 dx 3 = dx 1 dx 2 dx 3 (45) 

. However, as noted above there is no unique way to define the scale factor. In terms of the 
averaged quantity defined in (1441) we have 

a = \ (t pi +t P2 +t P3 ) (46) 
which, when we expand around t = 1, is approximately 

a(t«l)«i(2 + <) + 0(< 3 ) (47) 
. But in the limit of t small, we have 



a « l -t^ (48) 



, which is clearly a decreasing function; so the Kasner universe is not unambiguously expanding. 
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Redshift 



Redshift in a Kasner universe is given by 



(tR/tl 



{tR/tl 



\P2 







(t R /t E ) 



I'ii 



(49) 



1 

Z ~=3 



tn_ 1 
t E 



tn_ 1 
t E 



tR 

— ~ 1 

t E 



(50) 



In the circumstance when t R 3> t El 



. Of particular interest is the quantity AT/Tjj, the variation in CMB temperature from the 
average (accepting for the moment that the vacuum Kasner universe approximates a matter- 
filled one at a sufficiently young age), which is given approximately by 



AT 
Tr~ 



!• : I (I 1{P I. )''' 






3(t R /t E ) Pl - p3 -l 



-10 
-1 I eV 
2 





(t R /t E ) P2 

(t R /t E ) P3 


i{t R /t E r- pi -l o i 

2 



Vab 



(52) 



(except in the case when P2 = P3 = 2/3, in which event the (2,2) entry in ([52]) will read 2). The 
CMB in a mature Kasner universe has a pronounced anisotropy, with the observed temperature 
matching the average temperature only in a circle around the axis of anisotropy. Notably, the 
primary axis of the anisotropy is at a right angle to the axis along which the Kasner universe is 
contracting - not on a parallel axis! 



Hubble flow &; deceleration parameter 

The Kasner universe has Hubble flow 



H ab = 1 



Pi 











P2 











P3 



(53) 
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H = i^ 1 (54) 



and deceleration parameter 



Q 1 ;, = I ( o Ci -i>2 )/i>2 o i (■")) 





Q = -l (56) 

. In the limit that the parameter u^ooan observer in a Kasner universe would see a universe 
with a positive Hubble flow (redshift) over most of the sky, but see blueshift in a third direction. 
An observer looking only at averages, though, would not be able to distinguish between an 
isotropic universe and a Kasner universe merely by examining the Hubble flow; only with a 
complete picture of the sky is such a test possible. The Hubble flow in the case of minimal 
anisotropy has the form 



H ab {u = l) = - t \ 2/S | (.17) 



- appearing like a Friedmannian matter-dominated universe in two directions - and in the case 
of maximal isotropy 



lim H ah = \\ \ (58) 

u->0 t 1 



Similarly, an observer looking only at the averaged deceleration parameter sees a universe ac- 
celerating as though driven by a cosmological constant; only with good enough information will 
the observer notice a strong angular dependence in the acceleration field, which in the case of 
minimal anisotropy has the form 



-4 

Q b a {u = 1) = | 1/2 J. (59) 
1/2 

- decelerating like a Friedmann cosmology in two directions - and in the case of maximal 
anisotropy has the form 



lim Q b a = oo (60) 





10 GRAVITATIONAL WAVE NATURE OF BIANCHI IX 



25 



. Moreover, even though acceleration along two axes is negative in the least-anisotropic Kasner 
universe, the impact of the positive-accelerating direction is such that the isepitacn 10 ! of zero 
acceleration, the boundary an observer sees on the sky between regions where objects accelerate 
and objects decelerate, is a circle 83° from the axis of acceleration; only less than 8% of the sky 
appears close to "normal" to an observer expecting to record a Friedmann universe! 

While the vacuum Kasner universe is ruled out as a possible cosmology both for reasons of the 
CMB, which appears isotropic to a high degree[16], and due to the Hubble flow, which appears 
isotropic to the limit of the peculiar motions of galaxies below z = 0.3 [77], the surprising 
difficulties in distinguishing between its dynamics and that of a Friedmann universe serve as a 
reminder that sampling of cosmological parameters must be done in an unbiased fashion and that 
isotropy must be tested rather than assumed. The Kasner universe also has an application as 
a limiting case of the BKL universe[73], to which it appears identical for observers looking over 
a period of time that is small compared to the radius of curvature of the universe. Finally, the 
anisotropic Kasner universe serves as a limiting case to some types of cosmology described by 
the Bianchi IX model. 



10 Gravitational wave nature of Bianchi IX 



The Bianchi IX has been considered by cosmologists repeatedly since the establishment of general 
relativity to provide possible explanations for cosmological phenomena. 

Belinsky, Khalatnikov and Lifshitz discussed [74J a Bianchi IX cosmology (the "BKL cosmology") 
which undergoes several "bounces" as it evolves - rather than expanding from or converging to 
a point, it contracts along one axis while expanding along two others until the smallest metric 
component reaches a minimum value, at which point the axes swap roles. Misner[73 discussed 
a related form of Bianchi IX universe as the "mixmaster universe", pursuing an resolution to the 
horizon problem through the non-linearity of the Bianchi IX cosmology; through the mechanism 
of bounces, all parts of the universe are brought into causal connection. Bouncing vacuum cos- 
mologies are, like the vacuum Kasner universe, intrinsically highly anisotropic; while in the long 
run they tend to act isotropically due to the back-reaction of matter [751 176] they will still exhibit 
strong CMB anisotropyfljj]. Supernova data ([TJ [48] etc.) and CMB data on the value of CIm 
( \Sl\ etc.) coupled with the existence of high-redshift objects[72] rule out bouncing cosmologies 
to a high degree of confidence. 

The BKL cosmology undergoes anisotropic acceleration (see SECTION I9.2I) . Meanwhile, numerical 
modeling has suggested[98] that a matter-filled Bianchi IX universe will also undergo periods of 
acceleration. Therefore, we have good reason to suppose that a property of Bianchi IX may be to 
generate anisotropic acceleration, and that consequences of the Bianchi IX cosmology may reveal 
a dark energy candidate with none of the failings of scalar or exotic models. 

Wheeler showed|78j that an almost-isotropic Bianchi IX universe admitted a weak tensorial per- 
turbation that took the form of a wave (that is, solving an equation of the form f + nf (t) = g (t)). 
Grishchuk et. al. were able to generalize this result [TU]: 

The Bianchi IX space has frame vectors 

10 A neologism denoting a path of constant acceleration, similar to "isobar" or "isochor", from Greek "epita- 
chounse", acceleration. 
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e\ = (cos a; 3 , sin a; 1 sin a; 3 , 0) 

e 2 =(-sina; 3 , sin x 1 cos a; 3 , 0) (61) 



(i 



3 



: (0, cosa; 1 , 1) 



Consider the metric of a Bianchi IX cosmology: 



ds 2 = dt 2 - lah e1e)dx l dx3 (62) 

. We can split this metric up into an isotropic (Friedmannian) part and a non-Friedmannian 
part: 



ds 2 —dt 2 — a 2 F r\ ao e1 e^dx 1 dx^ — (7^ — a 2 F ri a b) e1^dx % dx^ — 

=dsl - (j ab - a 2 F r) ab ) e^e h j dx l dx j (63) 

. Grishchuk, Doroshkevich & Iudin showed that the object describing the space part of the 
anisotropic part of the metric at some moment in timeF^l 

G# ee 2 (ele) + efef ) - ^'W^ (64) 

, obeys the property 



= -(n 2 -3K)G$ (65) 

for n = 3 and K = 1; that is, Gfj is a tensor eigenfunction of the Laplace operator in a Bianchi 
IX space for waves with wavenumber n — 3. A similar property for open spaces is true of the 
Bianchi type VII models. p]5@ 

Lifshitz, in his development of the theory of cosmological perturbations [75l Ff9l \2\ ss. 115], claims 
that tensorial perturbations, including gravitational waves, can only have diminishing effect over 
time. Lifshitz is, however, considering only the class of local tensorial perturbations. 

In contrast, the gravitational waves in Bianchi IX will have wavelengths comparable to the radius 
of curvature of the universe. Kristian and Sachs note[25J that the wavelength of cosmic shear 
(and thus, all else being equal, of cosmological gravitational waves) must be at least 2 x 10 10 
years - longer than the Hubble radius[16j - and could potentially be far longer (see SECTION [TT]) . 

We will consider first the regime of weak gravitational waves, in an almost-isotropic universe, 
and then "quasi-isotropic" waves: that is, the regime in which components of the metric evolve 
at equal powers of t. 

11 ap has been scaled here to equal 1 

12 We could also choose to interpret Bianchi I as the degenerate case of a flat universe containing gravitational 
waves of infinite wavelength with n = 0. The Kasner universe, however, is not such a universe: all the anisotropy 
is governed by a single parameter, u, so the system has an insufficient number of degrees of freedom. The Kasner 
universe is more like the Taub universe [108]. 
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10.1 Einstein equations in the tetrad formalism 

For a metric g a p, let the space-space part of the metric be decomposed as in (|2"T|) . Similarly, the 
tensors 



Rij =R ab e a l e b J (66) 
Ta =T ab e^ (67) 

with all space dependence in the frame vectors. Therefore the Einstein equations can be rewritten: 



Rqo =kT 00 - ^kTg 00 

Roi =kT 0i - \kTg 0l (69) 



2 

2 



Rat =k ( T ab - ]-T lab ) (70) 



If we have energy-momentum tensor 



T l/tv = (jp + e)u fl ,Uu-pgv,u (71) 
T =e - 3p (72) 



then 



Too = (p + e) u u - pg 0Q (73) 
T i = (p + e) u Ui - pg i (74) 
Tab = (p + e) u a u b - pj ab (75) 

If we then choose a synchronous Gaussian reference system, as we always have freedom to do, 



.9oo =1 (76) 
90i =0 (77) 



so the Einstein equations read 



-Roo =k (p + e) uoUtk (p + e) u Q u - kpg 0Q - -k (e - 3p) (78) 
Roi =k (p + e) u Q Ui (79) 
Rab =k (p + e) u a u b - kpj ab - ifc (e - 3p) 7 ah (80) 
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. If we then demand that our coordinate system be co-moving with matter, 

u°=l 
u l =0 

then 



•Roo =^fc (e + 3p) 
Roi =0 

Rab =2 k (P~ e )lo,b 



Let 



and 



_ ld_ i 3 ; _ ld_ 

dab - 2 at 9«e a e 6 - 2 ^ 7af) 



d = d abl ab 

The Christoffel symbols associated with our metric then become (2j ss. 97] 



1 00 — 1 Oi ~ 1 00 ~ u 

r°- -d 

pi _ TJ 

1 0j -a^ 

pi -pi 

1 jfe — 1 jfc 



where f are the Christoffel symbols associated with the three-dimensional metric tensor 
The Ricci tensor can then be written as[51 ss. 97]: 



i?no — — d — d^du 



l00 — — " — a a a b 

Roi =0 

Rab =d a b + dd a b — 2d ac d b — P a b 



or explicitly [10] 



d + d b a d b a =-\k{e + Zp) 
d a b + ddab - ZdacdZ - Pab =^k(e-p) 7 a6 

d*C& =0 
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where is the three-dimensional Ricci tensor constructed from F* fc . 



10.2 The curvature tensor for Bianchi IX 



Grishchuk explicitly gives the curvature tensors for all Bianchi types, and a general method for 
easily deriving them, in [7J. These tensors can be stated in removable and non-removable parts, 
the removable parts corresponding to time-dependent rotations of the space. Let the symbol 



be 



(98) 



Then where 



r °ab — Tp^ {labd + Idab ~ Jbda) 



(99) 



(these are analogous to the Christoffel symbols of the full space, but with different symmetry 
properties) the non-removable part of the curvature tensor is given by 



L a b — 2r c a j 6 c j + 2r^ b r < | £l | c ] + 2r c ad r^ 



[be] 



(100) 



where square brackets around the indices indicate the antisymmetric part of the tensor; and the 
removable part is given by 



-v c C 



c u ba 



+ ^ (fa v b ~ fbVa) 



(101) 



and finally the curvature tensor 



Pab — L a b — bbeda ~ b ac d^ — bb a d 



(102) 



. In the co-moving case that v a = we can simply state P a b = H a b- For the particular case of 
Bianchi IX (the frame vectors (|6Tj) ) and the curvature tensor when v a = reads, for diagonal 
components: 



P" 



{lfgV f9 ) fg 
2 Il7cd|| 



be 
7 Vac 



7aflghV f9 V bh 

Il7c.ll 



(103) 



and for non-diagonal components 



P b a = -2 7 C V 



bedf 



hdf\\ 



lacldfT] v 



(104) 
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where ||7ab|| is defined as the determinant of "/ a i,. The Einstein equations show that when v a = 
the non-diagonal components of 7 Q b must be zero, so as a consequence of our Gaussian choice of 
coordinate system therefore we can without loss of generality write the metric for Bianchi IX 

7ii = a%e 2a 

722 = a 2 F e 2 ? (105) 
733 = a 2 F e 2 ^ 

with all other space-space components zero, so explicitly the the curvature tensor P a b for Bianchi 
IX reads 



P n =-e- 2& (-e 4Q + (e 2 ^ - e 2 ^) 2 ) e 2a (106) 

P22 =\e- 2S (~e^ + (e 2 ^ - e 2a f) e 2 ? (107) 

P33 =\er 2S (-e^ + (e 2a - e 2 ?) 2 ) e 2 ^ (108) 

Pab -0, a ^ b (109) 

and the contracted curvature scalar 

P abl ab = 2a F 2 e~ 2S [e ia + + - 2 (e 2a + 2 ^ + e 2 ^ + e 2Q+2 ^)] (110) 



11 Einstein equations for Bianchi IX 
11.1 Exact equations 

Let the symbol 5 = a + j3 + 7 for convenience as in Q37p. For our chosen metric, we have the 
auxiliary quantities 



d u = (ad + a 2 a) e 2a (111) 

d n = (d 2 + ad + 4aaa + a 2 a + 2a 2 a 2 ) e 2a (112) 

d\=H + a (113) 

d=3H + S (114) 



and cyclic permutations in a, /3,7 thereof for 22- and 33-quantities. The full Einstein equations 
for Bianchi IX read 
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+a F e 



2a+2f3 p 2a+2 7 j_ ,,2/3+2^ _ 
25 



-2(5 



2 4 ( 

_4« 



3 4/3 _ g 4 7 _ 3e . 



^ + 2- 

+2a 



-e 

2 



-2-25 

7 1 C 



+ a + 2£ 3a + 

4« _ (g2/3 _ p 2 7 ^ 



h d6+ 

o2<5 



+2a- p 2 e- 2S 



0>F _i_ 2^ 

+2ap 2 e- 25 



i+^+ff (3/3 + 

e 4^_ ( e 2 7 _ e 2a) 2 _ e 25 

4+7+^(37 + ^)+^+ 



; 4 7 _ ( e 2a_ e 2/3-)2_ e25 



fee 



3* M") 



. We can also define quantities as 
tensor: 



(115) 
(116) 
(117) 
(118) 



components of a gravitational effective energy-momentum 



ke„ 



- 2 ( e 2«+2/3 +e 2a+2 7+e 2^+2 7 ^ 



I + /3<5+ 
'- (e^' - e 2a ) 2 - * 25 



f (37 + 5) + 7<5+ 

e 4 7 _ ( e 2a_ e 2/3-)2_ e25 

2a + 6f^d + 2a 2 + a/3 + a A f - /3 7 + 

5e 2(a-//- 7 ) _ 3e 2(0-a- 7 ) _ 3e 2( 7 -a- 

+6e- 2a - 2e- 2 T - 2e~ 2/3 - 1 
3ff /3 + 2/3 2 + d/3 - d 7 + /3 7 + 

8-a- 7 ) _ 3e 2( 7 -/3-a) _ 3e 2(a-/3- 

V +6e- 2 ' 3 - 2e- 2Q - 2e- 2 ^ - 1 

2 7 + 6 ff 7 + 2 7 2 - a/3 + + £7+ 
_ 2 / 5e 2 ^^-«) - 3e 2 < Q -^-^) - 3e 2 ^^- 
^ +6e- 2 ^ - 2e- 2 ^ - 2e- 2tt - 1 



2/3 + 6^/3 + 2/3 : 

+ a F 



(119) 
(120) 
(121) 
(122) 
(123) 
(124) 
(125) 



(all of which are zero when a = /3 = 7 = 0). The Bianchi identity IT demands = pg = p' g 
so define the averaged gravitational pressure 



(3) 
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H 2) 



(126) 



25 



6^5 



2(d 2 + /3 2 +7 2 ) + 
+2 {af3 + cVy + £7) + 

_ e 2(a-f3-~f) _ e 2(/3-a- 7 ) _ g 2( 7 -a-/3) 

+2e~ 2a + 2e- 2 ^ + 2e" 2/3 - 3 



. Finally, 

fc (e g + 3 Pg ) = 25 + 4—5 + 2 f d 2 + /3 2 + 7 2 ) (127) 

Define a quasi-conformal coordinate 77 by ccft = apdrj; note that this fixes the relationship 
between t and i] up to the level of the characteristic length ai and a constant which can be set to 
zero. Given the impossibility of selecting a unique and objective definition for the scale factor, we 
do not define the conformal time using such a function. Define a correction term to the matter 
energy density q such that 

e = e F (l + q) (128) 

. In f^-time, the Einstein equations for Bianchi IX, subtracting background terms on both sides, 
read: 



-2(5 



a'/3' + a'Y + p'i + 2% S'+ 
2 ( e 2a+213 + e 2a+21 + e 2 ^ +27 ) - 



3 4 Q _ g 4/3 _ e 4 7 _ 3i 



,25 



a" + -£- (2a' + 5') + a'6' 



-2e 



-2<5 



,4a 



(e 2 ? - e 2 ~< 



) 2 -e 2S 



(3" + 2£ (2/3' + 5') + &8'+ 



+2e 



-2<5 



=4/3 



3 2 7 _ e 2a 



\ 2 -e 25 



1 " + ^(2 1 ' + 5') + i5'+ 
2e- 2S - (e 2a 



..2(5 



a F kepq 

-apkepq 
-apkepq 
-apkepq 



1 — w 
~ 2 

1 — w 



l-w o 



(129) 
(130) 
(131) 
(132) 



We also note the Einstein equations have an exact formal solution 



ke=(Sa'/e- s ) 1+w 



(133) 



where S is a constant of proportionality such that S 1+w has dimensionality of length to the 1 + 3w 
power. Finally, if we define the gravitational equation of state w g = p g /e g we see that necessarily 
w g = w and the Einstein equations can be read as 

13 Equation II127H corrects an error of sign in 1101 equation (27)]. 
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kPg + wapkepq = kp g 2 ' + wapkepq = kp^p + wapkepq = apkepq + ke g = (134) 

. In other words, the effective energy-momentum tensor created by cosmological gravitational 
waves equals minus the back-reaction on matter energy density and pressure. Note that the 
quantity ke g /q is necessarily negative. 



11.2 Solutions to the Einstein equations at zero order 

For convenience, define the variable x = 1+ ^ w rj. Then at zero order the Einstein equations 
for a Bianchi IX universe have, for arbitrary constant equation of state, the following solution 
and auxiliary quantities, which are identical to the solutions to the Einstein equations in the 
unperturbed closed Friedmann cosmology: 



ap =a,i (sinx) 1+3 "' 

1-3-w 

a F =a,i (sinx) 1+3 ™ cos a; 



l + 3w 



1 - 3k; 



l + 3w 



(sin x) 1+3 ™ cos 2 x — (sin x) 1+3n 



o! F jap = cot x 
Hp 



a,; 1 cot x esc x 



Qf = 



l + 3w 



sec 2 x 



(135) 
(136) 

(137) 

(138) 
(139) 

(140) 



. The quantity a; represents a characteristic scale for the universe, and in the background case 
represents the radius of curvature of the universe at the extent of its maximum expansion. We 
treat di as an arbitrary constant for the time being. 



11.3 Solutions at linear order 

To first order the Einstein equations take the form: 



2 — S 1 — 2Si =a F kepqi 
ap 



a{ + (2a[ +5' 1 )+8a 1 - 4Si 
ap 

aF 

7 f + ^(2 7 ( + £i) + 87i-4<5i 
ap 



1 




w 




2 




1 




w 




2 




1 




w 



a F ke F qi 
apkepqi 
apke F q 1 



(141) 
(142) 
(143) 
(144) 



where the subscript 1 denotes a first-order small quantity, that is, a quantity small such that in 
the first approximation its square is negligible. The formal solution f| 133[) gives us, to first order, 
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a%ke F qi = - (1 + w) S 1+u 'a^ 1 ^ 'Si 



(145) 



. Meanwhile, we can always choose to let S take on its Friedmannian value[TD], so S 1+w = 3 
Therefore: 



«; 



*1 + I — a l 
dp 

nil , n^F al 



dp 



II i n a F I 

h + 2 — 7l 

dp 



+ [3 (1 + w) esc 2 x-2]Si =0 

dp 



ap 
i 

Fx' 



a 



dp 



a 



dp 



1 


o 

— W 


to 




2 




1 


o 

— W 


csc 2 




2 




1 


O 

— W 







(146) 
(147) 
(148) 
(149) 



which gives us the solution: 



Si 



3+3ti 

c\ cos x (esc x) 1+3u 



(150) 



. The term governed by ci is a "removable" perturbation, that is, one not arising from a physical 
phenomenon but from small changes in our selection of the scale factor. Grishchuk, Doroshke- 
vich & Iudin argue[10j, and Grishchuk later proves in the case of high-frequency gravitational 
waves[103j, that the the removable perturbation arises in the transformation from i-time to con- 
formal time, and represents a small change in the value of 77 - that is, from the constant that 
has been implicitly set to zero in the relationship dt — dpdr] so that t = coincides with the 
metric singularity from a F (0) = 0. This coincides with the argument made by Bardeen |105| 
with regard to scalar and vector perturbations with wavelengths longer than the Hubble radius; 
Bardeen recommends a gauge choice minimizing shear. We always have the freedom to set c\ to 
zero but do not do so yet. In a radiation-dominated universe, we have 



SI 



radiation 



c radiation COSJ7CSc2j7 



(151) 



and in a matter-dominated universe 



SI 



matter 



c 



matter 



cos — esc 
2 



(152) 



Therefore the full first-order functions can be written: 



a" + 2 cot xa[ + 8ai = 3ci 



(csc X) 



3w 



3 + 3m 

(cscx) 1+3 ™ COS X 



(153) 



etc. Note that the right hand side contains no physical variables - no characteristic length or 
energy density. The Einstein equations at first order have the solutions (denoted with a tilde for 
the ci = case) 



1 1 EINSTEIN EQ UATIONS FOR BIANCHI IX 



35 



-radiation = ain3j] + c ^ cos3r/ ) CSC7] 

^radiation = sin3?? + c ^ CQs3 ^ cgC7? 

-radiation = sin3?7 + cos3?7 ) cgc?7 



and similarly for ,5,7 in a radiation-dominated universe, and 

-matter^ C<xi,i d sin3r; | C a2 ,i j cos 37? 

sin 77/2 drj sin 77/2 sin 77/ 2 di] sin 77/2 

etc. in a matter-dominated universe, in both cases constrained by the condition C a i t i + Cp\ t i + 
£71,1 = Ca2,i + C/32,i + C 7 2,i = 0. A general solution for any constant equation of state, in 
terms of orthogonal polynomials in a, exists but is far too cumbersome to be of practical use 
in this work. We introduce the notation C a \.\ etc. to be read in the following way; Cq2,i is 
an arbitrary constant associated with the function a, the hrst index denoting the mode of the 
solution (1 for growing, 2 for decaying), the second index denoting the order of the constant in 
an expansion assuming a, ft, 7 <C 1. For convenience, we will sometimes write a generic solution 
to the differential equation (|153| as 



&i = C Q ia?/i + C a2 sV2 (156) 
These solutions can be written in a less symmetric but easier-to-manipulate form: 



-radiation = c (2 cos 277 + 1) + C a2 ,i cos 3t? esc 77 (157) 



- matter 
ce-i 



C qM (16cos2?7 + lOcos?/ + 9) + 
-\C a i,\ esc 3 ^ (5 cos |?7 — 7 cos §77) 



(158) 



etc. When S — we recognize the homogeneous first-order Einstein equations as describing weak 
gravitational waves with wavenumber n — 3 and a wave equation of the form 



v" + 2 cot [x) v' + (n 2 - 1) v = (159) 

, in line with |10|'s description. In a radiation-dominated universe we have explicitly for the full 
first-order solution: 



radiation n sin3T 7 , n cos3t 7 , c i „„„ 2 / 1RfV , 

a 1 =C a i i- hC a2 ,i- h — cos 7/ esc r] (160) 

sin 77 sin t] 3 

etc. and in a matter-dominated universe we have 



matter _ d sin 3?? C a2 ,i d cos 377 d 17 3 77 

sin 77/2 d?7 sin 77/2 sin 77/2 0(77 sin 77/2 3 2 2 
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. It is common to refer to the decaying "cos" mode of these gravitational waves as "singularity- 
destroying" [10j, in that they diverge as 77 — > 0, which could seem at first to imply lim 7 a b — > 00. 

It is worth remembering that as the functions a, /3, 7 appear in the metric as exponents, that is, 
7n = cipe 2a etc., decaying functions are not necessarily "singularity-destroying" for the following 
reasons: 



• their divergence must overcome the convergence of the Friedmannian term, which in the 
case of weak waves will occur when w < 2/3 but not generally; 

• functions of the form e~ x v for x < 0, y < are non-analytic near x = 0, that is, they are 
not described by convergent Taylor series in that region. 

As C a 2,i + C/32,1 + C-y2,i = 0, either one or two decaying terms preserve the t = singularity 
when the removable perturbation is removed, in a manner analogous to that found in the Kasner 
universe, in the case of weak gravitational waves (although the price of this is a divergence later). 

When discussing high-frequency, localized waves, it is easy to define an amplitude of the waves 
by (for example) normalizing a root-mean-square value over the wave's period. In the case of 
cosmological gravitational waves however this procedure is not possible in an absolute sense due 
to the diverging character of the decaying mode. Fortunately, mathematical conditions on the 
relation of linear-order terms to quadratic-order terms revealed at quadratic order (see SECTION 
lll.4| cause the term "weak" to give itself an objective meaning. If we wish to normalize the 
growing modes, they have the following RMS values: 



..RMS _ 
V\ = 



(1 + 3w)tt 



(l+3u))ir/2 



yfdr) 



1/2 



(162) 



radiation, RMS n: „„„ N 
y x =V3 (163) 

„J natter ' RMS =V259«16.1 (164) 



It is interesting to note that in matter, the decaying "cos" mode of ai,/?i,7i has the same 
^-dependence as the removable perturbation; a cosmologist attempting to remove what they 
assume, based on an incomplete picture of the sky, to be a removable perturbation may inadver- 
tently be suppressing evidence of a gravitational wave! 

Finally, the gravitational energy-momentum tensor's (entirely removable) components read, to 
linear order: 



fee 



kp 



(i) 
flU) 



kp 



(2) 

3(1) 



kp 



9(1) 

(3) 
ff(l) 



=3 (1 + w) — o cos x (esc x) 1+3uJ 
-2>w (1 + in) —x cos x (esc a;) 1+3 " 



(165) 
(166) 



while the back-reaction of the gravitational waves at linear order gives us matter EMT compo- 
nents which vary from background by: 



<7i 



-3(1 + w) ci cos x (esc x) 1+311 



(167) 
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; when removable perturbations have been removed, first-order weak gravitational waves have no 
effect on the distribution of matter. 



11.4 Solutions at quadratic order 



The Einstein equations to quadratic order read: 



2 cot x8 2 +[3(1 + w) esc 2 x-2]S 2 



a'2 + cot x (2a' 2 + 5' 2 )+&a 2 - 4<5 2 



/?2 + cot x (2a 2 +S' 2 ) + 8(3 2 - 4S 2 = 



i 2 ' + cot x {2j' 2 + <5 2 ) + 8 72 - 46 2 = 



3csc 2 a; (i^L_2 



2 °1 



-I[jf-(af + /3f + 7 f)] + 

3±^csc 2 a; (- (1 + w) 5 2 

-ai«i + 8 G8i - 7i) a - 
-16a 2 + 16ai(5i - 4<5 2 

3 i=h csc 2 x (- (1 + iu) 5 2 + ii^5 2 



-^ + 8( 7l - ai r- 
-IQPl + 16/Mi - 4<5 2 

3 1=22 csc 2 x (- (1 + w) S 2 + ii^J 2 

-7^' 1 +8(a 1 -/3 1 ) 2 - 
-16 7l 2 + 167i<5i - A5\ 



(168) 



(169) 



(170) 



(171) 



. Taking the ^T equation QUI) first, 



2 cot x5 2 + [3 (1 + w) esc 2 x-2]S 2 = 



1 /^,'2 



2 \"1 



5? 



3°1 



/3f 

+4 (o?+)8?+^ 



(172) 



. The homogeneous part ($2 of course has the same form as at first order, representing a removable 
perturbation, so 



S 2 = C2 cos x (esc x) 



(173) 



The complete solution in integral form is 



1 3 + 3n 

8 2 = - cos x (esc x) 1+3 " 



3 esc x 



x < 



(l+w) 2 
2 

_IA'2_ 



-±(af+/3f 



/3? 



Hi 2 ' 



sec 2 x (sin a;) 1+3 ™ dr] + c 2 > 



(174) 
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Define the following pseudo- vectors and their Euclidean dot products 

(Cal(l) > £7/31(1), C 7 i(l)) = a 
(Ca2(l)> C/32(l), C 7 2(l)) =T 

(Cq1(1) + + ^7l(l)) = CT " £T 

(^a2,l + ^/32,1 + ^72,1) 
(Cq1,iCq,2,1 + C^i^C/32,1 + <7 7 i,i(7y2,l) 



T • T 
(T • T 



= a 2 



(175) 
(176) 

(177) 

(178) 
(179) 



«i + 0? + 7i -^ 2 y? + t 2 2/ 2 + 2 <T • ryiy 2 

~/2 1 o/2 1 — /2 2 /2 , 2/2,f, ' ' 

a<l + /3i +7! =cr y l +t y 2 + 2er • -r?/ l2 / 2 



(180) 
(181) 



. Note that the solution 62 = is excluded except in the case of a vacuum; this further excludes 
the definition 

a 2 = (711722733) 1/3 (182) 
as a useful definition of the scale factor. More careful examination reveals that the definition 



a 2 = -f ab V ab 



(183) 



is incompatible with the second-order equation as well. When all removable perturbations are 
set to zero, 



^on-removable = cog x (cgc x) m „ 




o 2 {2yi + yp 



(2?/! 

- (4yiz/ 2 + 52/12/2) 



tan 2 x (sin a;) 1+3tu dr\ 



(184) 



n We will discuss solutions to this equation term-by-term, noting that these terms can be solved 
entirely from information we obtained at first order F^l 



Contributions from the removable perturbations 



Contributions from the removable perturbations at second order have the explicit forms: 



14 The Einstein equations for weak gravitational waves in a Bianchi IX universe have the elegant feature of being 
integrable in closed form, always reducible to functions form sin {nrf) csc fc (rj) and cos (nrj) csc k (rj). Theoreticians 
working in regimes of higher-frequency gravitational waves in a fiat background may find it felicitous to approx- 
imate a Euclidean universe as a closed one in order to avoid mathematical inconveniences associated with the 
function sine (t) ! 

5 Li and Schwarz[107] obtain a similar result for a flat universe, but apply their results to a different domain. 
The averaging scheme they propose is not an applicable approach for cosmological gravitational waves. The result 
is generally stated in [2] ss. 96]. 
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In a radiation-dominated universe: 

^removable = _ £l ( 4 sin 2 \ + tan 2 \ + cot 2 \ - 2) esc 2 r, + c 2 cot r, esc r? (185) 

Note that the terms deriving from the first-order removable perturbation diverge as O (??~ 4 ), 
while those from the second-order removable perturbation diverge more slowly, as O (i] 2 )- 

In a matter-dominated universe: 



^removable 

°2 



2 

= - ^- ( 3 esc 4 - + 8 esc 2 - - 10) esc 2 - + c 2 cot - esc 2 - (186) 
12 V 2 2/2 22 v; 



. Similarly, terms deriving from the first-order perturbation diverge as O (77 6 ) and so at small 
77 will dominate terms deriving from the second-order removable perturbation which diverges as 

0(rr 3 ). 



Contributions from the growing mode 

Contributions from the growing mode have the following forms: 

^growing = ^ cQtx (csca , )T ^ j tan 2 x Q y / 2 + 2y 2^ (sina ,)T^ dr] (i 8 7) 
. We can already discern that the sign on 

growing 

must be positive in a young universe. 



In a radiation-dominated universe: 



£growmg,radiation = o^.^.^ cot v CSCf] ^^ cos3l] + ^ cos 5r) + 2 sec v 
note the diverging contribution of O (r)~ 2 ) from growing modes. 



(188) 



In a matter-dominated universe: 



*-V +™ sin 77 -^sin2r,+ 



^ r0WmS ' matter = -matter -t? esc 2 ? ( +^ sinS, - 33 sinV ] (189) 



2 2 



+ ^sin5?7 + 900 tanf 



. In contrast to the radiation-dominated case, the growing mode's contribution does not diverge 
in a matter-dominated universe (the term in brackets equals + O (i] 5 ))- Approximating to 
lowest orders in i], 

.growing, matter 2 f^ A r 2 21641 A 

5 2 S ' ~ -matter ~ ) (l90) 
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Contributions from the decaying mode 

In a radiation-dominated universe In a radiation-dominated universe, the functions yi and 
2/2 have the property 

y\ + y\ = esc 2 v (191) 
while the functions y[ and y' 2 are similarly related by 

y'i + V'i = (8 sin 2 77 + I) esc 4 rj (192) 
. This simplifies calculations as we can readily say 

^decaying = ^ ^ ^ * j^lA _ ^growing (ig3) 

; in a universe old enough that the diverging terms are negligible, the decaying mode intrinsically 
decreases the scale factor in the same way that the growing mode intrinsically increases it. 

In a matter-dominated universe In a matter-dominated universe, 



y 2 +y 2 = csc 4 |(9+icot 2 |) (104) 



and 



y'i + y'i = ^ csc8 \ (-608 cos r\ + 140 cos 277 + 477) (195) 



so we can state 



/ 18/7 + 2450 tan 2- \ 2 
5 decaymg,matter = ^ ^ V ^ V I 15705 cot n_l _ Ingrowing (ig6) 

2 2 V -^«*3«*3 ) ° 

. It is interesting to note that, due to the growing mode contribution's much slower contribution 
to change in the scale factor, the impact of the decaying mode on the dynamics of a young 
universe can be many orders of magnitude greater than the impact of the growing mode even 
when the decaying mode is several orders of magnitude weaker than the growing mode. The ratio 
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so in a matter-dominated universe with 77 w 10 1 the decaying mode will have a greater impact 
on cosmic dynamics as long as > lO" 8 ^^. 



Contributions from the <x • r term 



The contributions are described by the equation 



r mixed 



<r • x cos x (esc x) 1+3,1 



4yi2/ 2 + ^y'lVi ) tan 2 a; (sin a;) 1+ ~ 3w drj 



(198) 



and have the following explicit forms: 



Radiation-dominated universe In a radiation-dominated universe, 



r mixed, radiation 16 . . . 

5 2 = — cr • x radiatiQn sin 77 cos 77 (3 cos 2ti + 2) 



(199) 



Matter-dominated universe In a matter-dominated universe, 



rmixed *.V 2 V ( -4 cos 77 - 24 cos 2 77 - cos 3?7- 

<ft B = <T ■ T COt - CSC - 15 , , ' _ ' 

1 2 2 \ - ^ cos 47/ + 5 cos hr\ 



(200) 



Gravitational waves at second order 



Turning now to the R b a equations (|169l 11701 II 71 p . to second order, the Einstein equations for 
e — p( a '-terms read: 



a'2 + 2 cot xa' 2 + 8«2 
-3 [f (1 + id) esc 2 x + 1] 5 2 



-3 + ^ (1 + 3w) tan 2 x+ 

+ ^(l + w) (3to-l) + 
f (1 + w) 2 - fto) csc 2 x 
-oi*! + (6/3 2 - 16/3 l7l + 6 7l 2 ) 
-18a 2 + 16ai«5i 



(201) 



etc. If we suppress all removable terms, as we must for any practical observation of second-order 
terms, and taking into account (I184[) . this further simplifies to 



a'i + 2 cot xa' 2 + 802 — 3 



w 



^ (1 + w) esc 2 x + 1 



-26a 2 + 
-14/3? + 14 7l 2 - 



i K 2 



P? +1?) 



(202) 
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. Recalling the form of the gravitational waves including the removable perturbation at first 
order, make the simple transformation Q12 + to arrive at the equations: 



a 2 + 2 cot xa 2 + 8a 2 = 40 



(a\ +(3 2 1+ 7l 2 ) 



(203) 



etc; we recognize that linear-order gravitational waves act as a driving force on the waves at 
quadratic order. The solution of this equation is straightforward but tedious and we arrive at 
the following solutions: 



In a radiation-dominated universe 



„ radiation 

a 2 



C a l,2 

+40 (i^ 



sin 3r; , cos Zr\ 

r Q 2,2 



sin 77 



-40 



Ca2,l) 



sin rj 

sin 3r} 1 cos 3r/ 

6 ' sin rj 

1 sin 3)7 j 5 

24 



-401 



36 sin 77 

1 COS32 + 

b ' sin 

1 sin 5?? 1 

16 sin 77 6 

/ 1 „ sin 3?7 
6 1 ' 



2C a i ) it7 a 2,i) 



36 sin 77 

277— tt) cos 377 — 2 sin 377 ln(2 sin 77) 
sin 77 

^ cot 77- 



sin 77 

1 cos 377 1 cos 577 

36 sin 77 32 sin 77 



(204) 



etc. with the second-order constants C a \ 2 etc. constrained such that 



Gal, 2 + C/31,2 + C-yl,2 — Cal,2 + + C-yl,2 — 

To lowest order in 77 the solution for a 2 reads 



(205) 



radiation 

a 2 



C al . 2 (3 - 4ry 2 ) + 20 (±a 2 - C 2 al x ) ( 



0-1 



175 

36 



or" \ „— 1 , 1 rnon-removable 

^al,lW2,l )V + 3O2 



(206) 



etc. For the pure decaying mode, the contribution from S 2 dominates, while for the pure growing 
mode and the mixed term the contributions from the homogeneous parts of a 2 dominate. 



In a matter-dominated universe For a matter-dominated universe, the gravitational wave 
equation to second order has the following solution^!: 



a 2 — 



C a l,2CSc2j- 



7J _d_ sin 3?? 

sin rj /2 



Ca.2,2 CSC 



-a 



growing decaying mixed 



a 



rj d cos 3rj 
2 drj sin r//2 

ijnon-r 



+ 

emovable 



(207) 



16 There is no "royal road" to the explicit statement of this function, which was derived by substitution and 
variation of parameters with the assistance of a computer algebra system. With foreknowledge of the form of the 
solution, the equation l|203ft can be solved through the method of undetermined coefficients; this requires solving 
a 21-dimensional linear system. ^ 203 ft may also admit a solution through the method of Fourier transforms, but 
only under torture. 
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growing s g/'l 2 _ C 2 i] 



70 A-in=0 9n COS 717]+ 

/ -1128960 cos ^ + 
17 ^ +806400 cos ^ 
. +E"=o^sin(^ 



(208) 



,9o 



= 32900,gi = 443310,.g 2 = 90230,33 = 354221, g 4 = 20195, g 5 = 248918, 

36 = -57025,37 = 68911,3 8 = -37880, 39 = 15440,3m = -22400 
h = 1166543, hi = -1664285, h 2 = 888216, h 3 = 990580, h 4 = -1262310, h s = 677390, 

he = -363895, h 7 = 197841, h s = -116900, h 9 = 66864, h w = -34304, h n = 8960 



,f™^ „ fa* - Cj Ifl ) (82630- (209) 



, 1/1 \ , / -ftanlE^o-?" 008 ™^ 

4 6Caymg 4 ^-CL, ^ + jl *.C06»*f I 

- 1 +ln(-2sin 2 §)Eto^cos«7 7 



Jo 



= -34020, jx = -17010, j 2 = 153090, j 3 = 22680, j 4 = -113400 

fc = 58329, fci = -514422, k 2 = 368937, k 3 = 675396, 
fc 4 = -678540, k 5 = 31500, fc 6 = 61250 

Z = -5670, h = 102060, Z 2 = -73710, l 3 = -136080, h = 113400 



mixed- 4 ( 1 _ r- ^ \ ( § E«=o m « cos »?~ \ f -m 

12 = 105 la" ' T " CaMCa2 'V 2 I - cot i ELo "n — , J (2U) 



777.0 = 2310, mi = -39270, m 2 = -9240, m 3 = 46200 

n = -936, 77i = 15693, 77 2 = 30204, 773 = -58700, 774 = -25200, n 5 = 42000 



a mixed w _ 32 n . T _ Cal lCa2 \ v -s £ 77, 



n=0 



etc. The statement of the solutions to the gravitational wave equations to quadratic order in the 
matter-dominated universe are original to this work; the radiation-dominated quadratic order 
wave equations were presented in |10) . Note that J\ l n = Y) m n = 0. 

Most interesting is the presence of ln-terms in (j!93|) and (|210p , which on the one hand indicate the 
appearance of the power-law behavior of metric coefficients which typify the Kasner universe and 
the BKL universe in its quasi-isotropic phase; but which on the other hand show the breakdown 



1 1 EINSTEIN EQ UATIONS FOR BIANCHI IX 



44 



of our approximation scheme and the limit of regular perturbation theory in solving the problem 
to hand; the Taylor expansion of the growing mode of indicates further that waves must be 
very weak (||cr|| = O (l0~ 4 )) for the approximation scheme to be rigorously valid. In any case, 
indications are that the growing mode of hypothetical cosmological gravitational waves should 
be very much stronger than the decaying mode (see SECTION [19]) ; we will not need to make use of 
the second-order solutions for the decaying mode and from here on will treat the decaying mode 
as being linear-order weak, that is, C 2 21 ~ C«2 l ~ C72.1 ~ Cq2,2 ~ C$2,2 ~ C 7 2,2 ~ r 2 ss 0. 



11.5 Strong growing waves in the quasi-isotropic regime 

[TUl part 3] begins the development of equations for a radiation-dominated universe describing 
strong gravitational waves in Bianchi IX. Similar equations in a matter-dominated universe are 
useful in considering observed acceleration, as AQ sw — 1. 

Consider the equations (|129H132"|) . Assume a solution of the form 



71=0 

00 

/?=E c ^ 2 " ( 212 ) 

n=0 

00 

7 = E c ^ 2 " 

n=0 



with the terms constants. It is convenient to define e 2c ° = A, e 2c ° = B, e 2c ° = G. In a 
matter-dominated universe, to lowest two orders the solutions read 



7 «cj + 



20 



20 



1 



ABG 
1 

ABG 
1 

ABG 



(5A 2 - W 2 - 3G 2 + 6BG - 2AB - 2AG) 
(5B 2 - 3G 2 - 3A 2 + 6AG - 2BG - 2AB) 
(5G 2 - 3A 2 - 3B 2 + 6AB - 2AG - 2BG) 



(213) 



where Cq , Cq , Cq are arbitrary; if we want to preserve the Friedmannian value of S then we need 

<S + <£ + cg = (214) 

|10) . We always have the freedom to set one of these to zero by a simple scaling of the metric; 
this preserves the two degrees of freedom for the gravitational wave. 

If we apply the condition (|214|) and set the parameter Cq = by scaling, then the strong growing- 
mode waves are described by 
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eg eM (215) 

4 = - <% (216) 
cl =0 



(217) 



1 

20 

' 2 "20 

•v 1 



= 7 ^(-5A 2 + 2A + 6-6A" 1 +3A~ 2 ) (218) 
(3A2-6A + 6 + 2A" 1 -5A- 2 ) (219) 



,, . (3A 2 + 2A-10 + 2A- 1 +3A- 2 ) (220) 

with the single parameter Cg determining the whole system. Note that setting Cq = does not 
imply 7' = 0. We can also qualitatively say that for any value of A, two of functions a, /?, 7 will 
be positive, as will S, unless A — 1 (the background case), in the regime that Ar\ is sufficiently 
small that A 3 rf is negligible. 

The functions (|212[) are linearly independent with y ma tter ^ i owes t order in 77 and therefore 
can be used together to describe a matter-dominated universe with arbitrarily strong growing 
gravitational waves and weak decaying gravitational waves up to order rj 2 , as long as the series 
(|212[) converge. 



11.6 Dynamics 

As in the Kasner universe (see SECTION 19. lj) , it is useful to generalize quantities pertaining to 
the expansion of space which are spherically symmetric in Friedmannian cosmology. 

In terms of our statement of the metric (|105|) , the generalized dynamical quantities for our space 
are 



e a 

a ab = a F I | (221) 

e 1 



a = ia F (e a + e 13 + e 7 ) (222) 





a F /a F +/3 J (223) 

ap/ap + 7 



H = °* + U (224) 
ap 6 

f 'dp / a F -\- 2Hpa+ 
+(i + a 2 

dt" " L {Hp + af 




1 d 1 1 V +ct + a 2 
Q\ = ±H u nic-S\ = -^ — : -g (225) 



1 1 EINSTEIN EQ UATIONS FOR BIANCHI IX 



46 



etc. 



/ a F /a F -\-2H F a.-\-a-\-a 2 



V 



d F /a F +2H F p+/3+l3' ! 

(H F +f>Y 
I if /nF+2gF7+7+7" 



(226) 



/ 



. Our goal in undertaking the arduous task of solving the Einstein equations has been to derive 
the impact of long-wavelength gravitational waves on cosmic dynamics, particularly acceleration. 
We are now in a position to begin to discuss this impact. 



Let each quantity in section (jll.6|) be expanded out into a background term plus corrections, 
such that for example 



n ~ „(°) -Ln^ ,„(2) 

dab « a ab + a ab + a ab 



(227) 



Then the zero-order, background terms are simply 

(0) 

a ab = a FVab 



H 



(0) 



'-H F r)ab 



(o)n b 



Q°a =QFS° a 



(228) 
(229) 
(230) 



. While the gravitational energy-momentum tensor vanishes at first order with the removal of 
removable perturbations, the presence of weak gravitational waves can affect observed dynamic 
quantities. At first order: 



H 



0(1 



(i) _ 














A 











71 








) = 












Oil 






















7i 



(231) 



(232) 



(233) 



(i) 



H(i) - 3*1 



Q\ = -Hp 1 [2 (Q F + 1) ax + Hp 1 ^} 



(234) 
(235) 



etc., 
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1 



2(Q F + 1)S 1 +H F 1 5 1 



(236) 



Thus we illustrate the need for truly representative sky coverage in considering the problem of 
acceleration: gravitational waves can contribute to anisotropic acceleration even when they do 
not affect the distribution of matter. In domains when the first derivatives of a wave is small 
(that is, near peaks and troughs of the wave), the accelerative effect will not be accompanied 
by a large change in the Hubble flow. As before, a failure to completely suppress the removable 
perturbation may lead to incorrect evaluation of the strength of decaying modes. To first order, 
non-zero contribution to the average over the whole sky of the perturbations is removable; first- 
order weak gravitational waves in Bianchi IX do not produce isotropic acceleration. 

To quadratic order, the dynamic quantities have the forms 



(2) 

a ab = a F 



a 2 + aj/2 

+ 

7 2 + 7l 2 /2 



(237) 



«2 = gOF 



5 2 + -(al+Pl+ 1 l) 



(238) 



rr(2) _ 



a 2 






















7 



(239) 



(2) 



(240) 



{2) Q\ ~ 



= -HI 



2(Q F + l)a 2 +H F 1 a 2 - 
-3H- 1 (Q F + 1) a\ - 2H~ 2 a l a l 



(241) 



etc., 



Q(2) 



-H, 



2(Q F + l)d 2 + H F 1 S: 
-2H~ 2 (diai + /3i/3i + 7171 



■7i 2 



(242) 



. At second order we begin to see a consequence of the non-linearity of the Bianchi IX Einstein 
equations which is potentially very important in the study of cosmic dynamics: isotropic changes 
to the Hubble parameter and to acceleration from anisotropic metric terms. With our knowledge 
of the Einstein equations at first and second order (|153ll72l202p we can show this explicitly: 



(2) 



(3w-f 
-40tana;o!x 



tan x < 



tan x 



+ 3w) tan 2 x) a' 2 — 8 tan xa 2 — 

(l — 3w — (1 + 3w) tan 2 x) t&nxa' 2 

+16 tan 2 xa[ai + 
3 [f (1 + iu) csc 2 x+l] S 2 + 



-14(a 2 +/J 2 + 7 2 )-iK+/?f +7f) 



(243) 
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etc. and 



>{2) 



1 2 

= — tan x 



\ (1 + 3w) 2 sec 2 x£ 2 — 
-2 (1 + Zw) sec 2 x (a\ + $\ + j 2 ) + 
-\ [l + 15to + 5 (1 + 3u>) tan 2 x] (of + ^ + 7 f ) 
— 16tanx (a^ai + fi'^i + 7(71) 



(244) 



. Isotropic acceleration with quadratic-order strength arises from the non-linear interaction of 
linear-order gravitational waves; but in the regime of \a\ , , | 7 | <C 1 the gravitational waves at 
linear order will dominate measurement of cosmological parameters. 

In a matter-dominated universe with 7/ small the deceleration terms become, defining 



AQ a b = Qt- Q F 5 a b 



(245) 



AQ = Q-Q F 



(246) 



AO 1 

^^1, matter ~ a 



tan 2 (a' x + a' 2 ) — 8 («i + a 2 ) 



-40af + faf + 16 tan §aia r 



+3<5 2 + 14 (a 2 + /? 2 + 7l 2 ) - I « + fi' 2 + 7?) 



(247) 



AQ 



matter 



48^ 



2<5 2 -8(a 2 + /3 2 +7i 2 ) + 
+ (a' 2 +fi' 2 +1 ' 2 )- 
-16 tan § (ai«i + /3^i +7W1) 



(248) 



etc. Explicitly, these will have the lowest-order forms: 



matter 



' 4 



C alA (280 - 259t7 2 ) + 
+C Q2 ,i (16tT 3 - 727?- 1 + ^firy) 
+C 2 al ! (710040 - 4687753ry 2 ) 4 
\ +a 2 '(-609590/3 +^|405r? 2 ) 



(249) 



AQ 



matter G 



^rj 2 (-4900 + 2983T7 2 ) 



(250) 



. These results are encouraging as if we choose ||cr|| ~ 10 4 (in order to make the gravitational 
waves weak) and rj ~ 10~ 2 to match (fT7j) . we obtain AQ} ma f+ er ~ — 10~ 6 , which has the right 
sign as well as all the contributions at both first and second orders going in the "right" direc- 
tion, toward acceleration. It is particularly encouraging that both growing and decaying modes 
contribute to acceleration to their lowest orders in 77. 



1 1 EINSTEIN EQ UATIONS FOR BIANCHI IX 



49 



11.7 Back-reaction 

Of interest in discussing the problem of acceleration is the effective equation of state of the 
gravitational waves' contribution to the energy density. Empirically, the equation of state of 
dark energy seems to be close to wx — — 1 (see SECTION [5]), where the quantity w x is related 
to the source of the energy such that the source evolves with regard to the scale factor at a 
rate of a~ 3 ( 1+Wx \ As noted in (section 19. ip there is no unique way to define the scale factor; 
but a condition of quasi-isotropy is that expansion in every direction in the current epoch is 
proportional, that is, that they evolve as the same power of time. If the decaying mode of 
the cosmological gravitational wave is weak, then this evolution will be proportional to the 
Friedmannian scale factor. 

To quadratic order, pi9p reads 

fce<, 2) = 3 (1 + w) a^ 2 esc 2 x5 2 (251) 

and so by (JM) 

9(2) = - (1 + w) a F 2 S 2 (252) 

. When the growing mode is dominant, 5 2 is always positive in a matter-dominated universe, and 
therefore q^ 2 ) is negative and thus the back-reaction appears to have negative energy density; a 
significant "mixed" cr • t term, however, can easily introduce intervals where q^) > 0. 

In a matter-dominated universe and when the growing mode is dominant, q^) V~ 2 which, if 
the universe is evolving with a scale as ap oc rj 2 . implies an equation of state for the back- reaction 
of wx = — 1/3 (as compared to an equation of state for a cosmological constant of wx = — !)• 
While no investigation of the equation of state of dark energy includes this value within its highest 
confidence interval, measurements of wx remain tentative, with large errors and high sensitivity 
both to single data points and to the algorithm for curve-fitting models to the data (see SECTION 

IS). 

The dominant term in (|119[) is the a'p/ap-teim. This stands in stark contrast to the commonly- 
considered case of gravitational waves in a background so slowly moving compared to the period 
of the waves that hp rs 0, in which instance the quadratic combination of first-derivative terms 
dominates. 

In regimes of stronger growing-mode gravitational waves, though, the scale factor as defined in 
(|39| will be more dominated by terms of higher, even order and so a a {, oc r/ 4 or higher. As the 
growing mode increases in strength, the equation of state decreases asymptotically toward a limit 
of wx = — 1 ; if the scale factor grows as rj 2s , the equation of state for the back-reaction is given 
by wx — (l/3s) — 1. As acceleration is empirically Qq = —0.6, this implies that in real life the 
gravitational wave strength is of order unity and therefore the effective equation of state is close 
to —1. Thus, the quasi-isotropic Bianchi IX model with strong growing- mode gravitational waves 
and weak or zero decaying-mode waves is compatible with the observed data on the equation of 
state of dark energy, without the invocation of a cosmological constant; the theory would be 
invalidated by definitive measurements of wx < — 1. 

In any case, the fact of wx < allows us to draw a conclusion regarding cosmic evolution. |76| 
notes Kasner-like cosmologies go through two stages of evolution: 
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1. A "vacuum" stage, where matter's influence is, due to its evolution as a -4 , weak compared 
to the influence of the anisotropic expansion and contraction, influence which, in light 
of (section [TP)) , we now understand to be the result of gravitational waves in the BKL 
universe; 

2. a "matter" stage, where expansion isotropizes[106j and is driven by, first relativistic (w = 
1/3), then cold, non-relativistic (w — 0) matter. Formally, the contribution of curvature 
to cosmic evolution becomes important in this era (wk = —1/3), but as the influence of 
curvature will be isotropic in Bianchi IX and the radius of curvature is very large compared 
to the Hubble radius (see SECTION fl7j) , curvature will not have a practical influence on 
observations in and of itselfl3 To this second stage we can add a third stage: 

3. A "dark energy" stage, in which growing modes of the cosmological gravitational waves 
which drove the initial isotropy return as the dominant influence on cosmic evolution. 



11.8 Amplification of gravitational waves 

Grishchuk observed[8j that when the background of a cosmology containing gravitational waves 
varies rapidly, weak gravitational waves can be amplified where they would otherwise, in a slowly- 
moving background, decay rapidly [79J. With regard to the Bianchi IX cosmology this is significant 
as when the growing mode of a cosmological gravitational wave dominates, the leading term 
in the gravitational energy density is of the form (a' F /ap) S' 2 = O (constant). Cosmological 
observations (see section [T7|) indicate the universe has rj < O(l0 _1 ). In this regime, the 
term a' F /ap = cot (77/2) 2/ry, which is dependent on the rate of change of the background, 
is arbitrarily large; therefore, weak waves may have an effect orders of magnitude greater than 
their amplitude. Similarly, the decaying mode of gravitational waves can have prominent or even 
dominant power in a sufficiently young universe even when the amplitude of the decaying mode 
is smaller than that of the growing mode. 



12 Conclusions 

Solutions have been presented for the gravitational wave equation for a Bianchi IX universe 
perturbed to quadratic order from the closed Friedmann case. Quadratic order is the limit of 
perturbation theory's applicability to explore nearly-Friedmannian Bianchi IX when decaying 
modes are sufficiently strong that they are not negligible. 

At quadratic order, the non-linear interaction of the gravitational waves produces isotropic 
changes to dynamic quantities. While this isotropic change is likely to be dominated in any 
particular direction by linear-order contributions from the gravitational waves, in the regime 
of strong gravitational waves they will become more important and potentially even dominant. 
Where |98| discussed the possibility of acceleration in a non-vacuum Bianchi IX universe only 
qualitatively, we have shown it explicitly as well as illustrating a clear link between acceleration 
and the gravitational waves which are intrinsic to Bianchi IX in its full generality. 

17 Formally we can also say that, due to the action of proton decay and positron annihilation, after sufficient 
time the w = phase will return to a w = 1/3 phase where the universe is filled with neutrinos and photons. 
Following this period there will be another return to w = as these free particles are absorbed by black holes. 
As these black holes evaporate by the process of Hawking radiation, there will then be a final return to w = 1/3. 
I80| gives a popular-science presentation of the universe in these phases, but as it was written only shortly after 
the discovery of acceleration its treatment of dark energy is highly speculative. 
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It is curious to note that the order- rj 2 approximation we have made in (section 111. 5|) . a and S 
in the normalization we have chosen take the form of Alexander polynomials^lO!}, 110J, although 
not Alexander polynomials for any knot of fewer than 11 crossings. Whether this mathematical 
observation is significant or coincidental is a subject for further debate, but as gravitational waves 
in Bianchi IX are moving equatorially around our background 3-sphere|10|. and as a sub-class of 
knots (the "torus knots") are constructed by wrapping one 2-torus around another it is conceivable 
there could be a connection. 

Back-reaction from growing modes of the gravitational waves appears to have negative energy 
density and an equation of state compatible with that observed for dark energy, especially in the 
regime of strong gravitational waves and quasi-isotropic expansion; when gravitational waves are 
strong, they become the dominant contributor to the evolution of the cosmos in an era following 
the era of matter domination. 

Therefore, from the perspective of cosmic dynamics, cosmological gravitational waves in a quasi- 
isotropic Bianchi IX universe are a viable candidate for dark energy, without the invocation 
of a cosmological constant and without requiring any modification of the theory of relativity. 
An analysis of the impact of these gravitational waves on the cosmic microwave background 
is necessary in order to determine whether constraints from the CMB are compatible with the 
observed data on acceleration. 



Part IV 

The Cosmic Microwave Background of a 
Bianchi IX universe 

While long-wavelength gravitational waves can cause acceleration in a Bianchi IX universe, the 
effect of such waves must be compatible with the observed cosmic microwave background in order 
to represent a practical model for explaining observed acceleration. 

Sachs & Wolfe initiated|23j the systematic study of the effect of perturbations on the CMB, 
following a formalism developed by Kristian & Sachs [25]. Sachs & Wolfe's work developed the 
theory of scalar, vector and tensor perturbations on the CMB in a flat almost-isotropic universe 
to first order. 

Sachs & Wolfe's work was generalized by Anile & Motta[2E| to the almost-isotropic closed and 
open Friedmann cosmologies, again at first order. While Anile & Motta begin to consider the 
impact of long- wavelength gravitational waves on the CMB, they choose to explore the impact of 
waves with scales much smaller than the Hubble radius. Anile & Motta subsequently[57] ruled 
out the existence of these waves at significant strengths in the observable universe. 

Doroshkevitch, Lukash & Novikov considered the impact of an anisotropic universe on the CMB 
in the case of the Bianchi VII, VIII and IX models [19], and concluded that a Bianchi IX model 
was potentially "compatible with observations, only if there was some secondary heating of the 
intergalactic gas". Doroshkevitch et aVs most important calculations are carried out on the 
assumption, then widespread, of VLm ~ 1 and as such are of limited applicability; interestingly, 
though, in their conclusions they note that if Hm < 1, "AT/T will be close to the maximum 
value only in a small 'spot' with an angular size 9 f=s 4f2" (where by "small" they give the example 
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of n M ~o.i => e»«23°). 

Sung & Coles analytically and computationally explore the impact of various unperturbed Bianchi 
models, including Bianchi IX, on the CMB[5T]. They report the useful theorem that "a gravi- 
tational field alone is not able to generate polarization", but do not consider the general case of 
Bianchi IX, only the isotropic case equivalent to the closed Friedmann universe. 



13 Geodesic equations 

The effect of the metric on the CMB is determined by examining the change in geodesies of light 
rays relative to an isotropic, background case. Let the subscript E denote a function evaluated at 
the time of the emission of a photon, and the subscript R denote that function evaluated at the 
time of the photon's reception. Then the change in the temperature of the background radiation 
T is given by 



. Consider the path of a light ray; let this be a four-vector denoted by fc M such that k^k^ = 0, 
with the light ray received in the direction k R — e l . The geodesic equation for the time part of 
fc M in a Bianchi cosmology reads 



Tr/Te 



1 



(253) 



z + 1 



- + rV^ = o 



(254) 



and the equations for the space part of the vector read 



dk a 

-tv + r oo + + T^k* + T a bc k b k c = 



(255) 



Recalling (|55J and fTUS]) the Christoffel symbols 




(256) 



The symbol e a i, c represents the Levi-Civita symbol defined such that £123 = 1 
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with all others zero; note that the form of the rotation coefficients guarantees that only anisotropic 
parts of the metric tensor will have an effect on k l (and therefore (5-terms, whether removable 
or non-removable always vanish in the geodesic equations; recall section 19. ip . Using the same 
method of conformally-related objects as described in ( |23| part He]), define the vector fc M : 
a 2 F k^ = fc M and the tensor j a b : a 2 F ^ a b — 7 a i>; recall that k R — —k R kf = 1. This gives us geodesic 
equations: 



dk° 1 



elk 



- + ( e 2 7 -2a _ e 2/3-2a) pp =Q 

dX v ' 

+ , 2a _ 2f} _ g27 _ 2 /n p p =Q 
dX v ' 

+ ( e 2^-2 7 _ e 2a-2 7 ) pp =Q 

dX 



(258) 
(259) 
(260) 
(261) 



. Despite the symmetry of these equations, their nonlinearity has inhibited the discovery of 
exact solutions and research into their properties is ongoing; see for example [24]. However, with 
solution up to quadratic order for the metric in hand (|160ll61|206|207j) . we can explicitly solve 
the equations in the case of weak waves. Let k a = k R + Ak a (A). Expanding out the geodesic 
equations to second order in the metric: 
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(262) 
(263) 
(264) 
(265) 
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. To first order, the equations are trivially solved by choosing A = r\ as the affine parameter; the 
problem of determining dX/drj is overcome by our choice of reference system, the lack of vector 
perturbations and the homogeneity of space: 
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. The relationship (|270[) explicitly shows the quadrupolar nature of changes to the CMB alluded 
to in [19J. An unremoved removable perturbation changes the temperature of the whole sky 
isotropically; this confirms the effect noted by Hwang & Noh|42|. 

The equations for quadratic-order corrections read 
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which due to the cancellation of the terms in the right column integrates trivially to 
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(reiterating the quadrupolar character of the change to the CMB, but generalizing it to anisotropic 
expansion); meanwhile for the space part of the vector 
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14 Redshift and CMB variations 

The geodesic of a light ray is related to its observed redshift by 

(*»<V)l5 

|23) . Having determined uq = 1 and ui — this simplifies to 
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. Meanwhile, the temperature field 
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15 Comparison with the observed CMB 



Five-year and seven-year results [16j [18] from WMAP[30] give the best picture to date of the CMB. 
The WMAP observations reconfirm the constraint of the quantity AT/T < 10~ 4 ^\; any change 
to the CMB from acceleration must be equal to or smaller than this value in order to be compatible 
with observations, placing an additional constraint on cosmological models. This implies that 
in the current epoch, and in the absensce of further special alignment, \a\ , |/J| , | — ^ | < 10 -5 . In 
a matter dominated universe, under ordinary circumstances, this implies (since rj < 10 _1 ; see 

SECTION [P7f 

< 1(T 5 => |C Q i,i| < 1(T 6 (287) 

< 1(T 5 =>> \C a2 ,i\ < 10" 8 (288) 



n „ .matter 



n „ .matter 



; meanwhile in a radiation-dominated universe, 



n „. radiation 

L'ai.iyi 
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W2,i2/2 ^ 1U 
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\c«iA < 10 



-5 
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(289) 
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. The coefficients associated with the decaying mode are constrained to be smaller than those 
associated with the growing mode without further theoretical considerations. 
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15.1 CMB anomalies 

Since the publication of the latest generation of CMB maps[28j, numerous claims have been 
made (for example, [20 OSES [39]) of anomalous structure in the CMB. While the WMAP team 
argue[r7] that these phenomena are not of statistical significance, if a quasi-isotropic Bianchi 
IX universe could produce any of the perceived patterns it would point the way toward further 
observational studies of the CMB to determine cosmological parameters, and establish the quasi- 
isotropic Bianchi IX universe as a viable model for cosmology. 

In all cases, we emphasize that the most likely explanation for any perceived pattern in the 
CMB which is not shown to be statistically significant is the null hypothesis: that is, the human 
perceptive phenomenon of pareidolia, the same phenomenon responsible for observing familiar 
shapes in clouds or the "Man in the Moon". 



Cold spots, "fingers" and the "Axis of Evil" 

Two compact, supposedly anomalous areas of low temperature have been noted in the CMB, the 
so called "cold spots". 

The first of these (called Cold Spot I in [17]) is a region[3T| covering approximately 15000 square 
degrees in the direction of the galactic center, much of which is 194 microkelvin|28j colder than 
the CMB mean temperature (AT/T R = -7.12 x 10~ 5 ). 

Particularly noteworthy regarding Cold Spot I is its membership in one of four "fingers" spaced 
at roughly 90-degree angles around the galactic equator, intersticed by four areas of higher 
(AT/T R = 7.12 x 10~ 5 ) temperature© Qualitatively, such a pattern is roughly consistent with 
the expected pattern if two of the functions a, j3, 7 > and if two of the the principle axes of 
the metric tensor lie on the axes of the cold and hot zones (implying the third axis points along 
the "Axis of Evil", see below). The so-called "Cold Spot II" reported by Vielva et al-ISH [3Z] also 
forms part of these "finger" structures [17J. 

Cold Spot I also has the angular size [19] predicts for the observed value of f2jvf ~ -3. 

Due to the coincidence of the cold spot with the direction of the galactic center, there are no 
optical observations in its direction (see FIGURE [S]), and therefore there is no data on cosmic 
acceleration in the direction of Cold Spot I. 

(Equation 1286(1 implies that any cold spot resulting from anisotropy in the metric should be 
accompanied by an identical cold spot at a point antipodal to the original spot. Tegmark's 
examination [55] of the one-year WMAP data on the CMB low-order multipoles revealed an 
alignment between the CMB quadrupole and octupole in the direction of (l,b) s» (—110°, 60°) 
along which the quadrupole is nearly zero, an axis which Land & Maguiejo found[32j extended to 
the 16-pole and 32-pole as well; the alignment has been dubbed the "Axis of Evil". While exam- 
ination of the three- year WMAP data[33j found the Axis of Evil to be of lower significance than 
initially thought (94%-98%), it still persists; the WMAP team's discussion of the alignment [171 
pt. 7] admits the "remarkability" of this alignment and, while assigning its existence to chance, 
does not attempt to explain the "Axis of Evil" in full. 

The Axis of Evil, which in equatorial coordinates [35, p. 43] lies close to RA 10:44 Dec +7.6°, 
falls within the zone in which redshift data has been collected for measurement of the cosmic 



The CMB dipole is defined as such a way as to be traceless, so J ATq Uac j rU p j e /TjidS = 0. 
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deceleration parameter. To simplest linear approximation with a pure growing mode, (that is, 
that the functions a and a' are both small such that a 2 ~ 0) this alignment rules out a CMB 
arising from cosmological gravitational waves as a source of cosmic acceleration. However, the 
fact of the alignment of the quadrupole, octopole, 16-pole and 32-pole indicates that non-linear 
contributions of gravitational waves to acceleration are not ruled out. 

The question of the overall magnitude of the quadrupole, which is only 14% of the expected 
value [281, l38j - has also been raised. The WMAP team|17[ pt. 4] agree with Tegmark that 
the depressed quadrupole falls within the 95% confidence interval for simulations of the CMB, 
but do not attempt an explanation for the unusually strong octopole term. Long-wavelength 
gravitational waves can easily explain both through judicious choice the arbitrary constants 
Cq.1,1 etc. in a manner compatible with the CMB. Efstathiou|29| supposes that the depressed 
quadrupole could be an indication of a closed universe; however the relationships he proposes 
generate zero contributions to the CMB power spectrum from the genuinely cosmological, intrinsic 
n = 3 waves found in Bianchi IX, and any observational test using his framework must rely on 
correct evaluation of gauge terms whose effective wavelengths must be far longer than the cosmic 
horizon. Furthermore, Efstathiou's conclusion that a closed universe would automatically require 
a scrapping of current inflationary models is contradicted by others; for example Guth argues 
that a universe that is closed but with a very large radius of curvature is not ruled out|40). 

The quasi-isotropic Bianchi IX model cannot provide an explanation for hemispherical dipole 
asymmetry claimed by Ericksen et. aZ.|39|. 



16 Conclusions 

The long-wavelength gravitational waves intrinsic to a quasi-isotropic Bianchi IX will cause a 
change in the cosmic microwave background with a distinctive quadrupolar signature. A radially- 
symmetric pattern of light deflections in the CMB resulting from shear may also be observed. 

The almost-isotropic Bianchi IX model can be compatible with the CMB as observed, and can 
provide an explanation for perceived anomalies observed in the CMB by COBE and WMAP. 
However, the existence of these anomalies beyond the level of statistical noise is not certain; a 
possible route of cross-disciplinary research is open in the form of examination of the phenomenon 
of pareidolia as applied to the CMB. 

Models of quasi-isotropic Bianchi IX relying on pure growing modes or pure decaying modes 
of the gravitational waves cannot simultaneously explain observed cosmic acceleration and the 
observed cosmic microwave background. Research into the non-linear regime of the Bianchi IX 
cosmology may elucidate the existence of a model of an accelerating universe in which isotropy 
is almost preserved. 
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Part V 



An accelerating Bianchi IX universe 
preserving an almost-isotropic CMB 



In order for a Bianchi IX universe to both appear nearly isotropic in the cosmic microwave 
background and to accelerate through the existence of long-wavelength gravitational waves, it 
must fulfill two conditions. The first is that the function k° (rjn) must have absolute value less 
than the limit imposed by observations of the cosmic microwave background, AT/Tr. The 
second is that at least one of the functions Q b a < 0. It is possible for both these conditions 
to be simultaneously filled while remaining compatible with other observational constraints on 
cosmological parameters. 

The idea of long-wavelength gravitational waves causing anisotropy in the CMB has been pro- 
posed, but not applied to the Bianchi IX universe. Grishchuk & Zel'dovich consider the possi- 
bility of long-wavelength gravitational waves existing in a Friedmann universe without violating 
the limits imposed by the CMB[41J, but do not apply their work to the gravitational waves 
of cosmological character which appear in some homogeneous cosmologies. Campanelli et. al. 
suggest that such a universe could exist and propose a Taub-type Bianchi I universe which also 
includes anisotropic dark energy as an initial explanation for the observed CMB, complement- 
ing Rodrigues[113j. Critically, they do not consider gravitational waves as a generator of the 
anisotropy and treat the parameters of the Taub universe as if dark energy were simply estab- 
lished by fiat. Similarly, Kovisto and Mota [TT5] do not look beyond the Bianchi I model and 
instead fall back on exotic theories to explain dark energy. 



17 Cosmological parameters 

WMAP[18] [16] has produced an all-sky survey of the CMB which, if the universe is almost 
Friedmannian, can be used to constrain cosmological parameters. 

Let the radius of curvature ao and conformal time r\ of the background Friedmann cosmology be 
treated as a free parameters; assume a closed universe. The WMAP seven-year data gives 



(WMAP's analysis includes the value of VIk measured by baryon acoustic oscillations reported 
in |101] V The radius of curvature, Hubble parameter and curvature energy density are related 



H =70.41 J-jkm/s/Mpc 
Sl K =- .0025 ±0.0109 



(291) 
(292) 



by 




(293) 



while the Hubble parameter, radius of curvature and 77-time are related by 
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H a = cot (no/2) (294) 

. Therefore we have limiting values (as defined by the 95% confidence boundary of the WMAP 
observations) 



a >1.12 x 10 29 cm (295) 

ry <0.0266 (296) 

and highest-confidence values 

a =2.68 x 10 29 cm (297) 

ry =0.00499 (298) 

. Meanwhile, the ratio of Hubble radius to radius of curvature is at least 

H a a > 8.67 (299) 

with a best-fit value of 



H a = 20.0 (300) 

. In other words, if the universe is closed, then the cosmological gravitational waves of the Bianchi 
IX cosmology are of much, much longer wavelength than the observable universe. 



Finally, from the value of the redshift of decoupling, z^g^ scattering = 1090, we can say by ([19 
that 



m/VE « 33.0 (301) 

. As the available data, including that from supernovae (see table [TJ, does not exclude a flat 
universe, we are always free, in developing the theory of Bianchi IX and acceleration, to set the 
parameter rj as close to zero as necessary; doing so will not, in and of itself, violate observations, 
but will instead be constrained by the impact of the decaying mode of the gravitational waves 
on the CMB. 



18 Compatibility with the redshift 

Of all the observed cosmological parameters observed by WMAP and other probes of the CMB, 
the ones that are directly observed are AT/Tr and ^ as ^. sca ttering' F rom these we can say that 
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in the current epoch the universe appears isotropic and that its expansion since last scattering 
has, on average to the present time, been isotropic; neither of these facts necessarily imply that 
the overall expansion was isotropic at any time before the present. Instead, the condition of 
quasi-isotropy simply implies that 



+ ^ ab ,ok a k b « 
or] 2 



(302) 



This implies that shear is small, so 



k a «A)q 



z + 1 «of (jin) /a F (rj E ) 



(303) 
(304) 



as in the background Friedmann case. 

We can obtain a near-zero value to the wave functions in the present epoch by admitting the 
presence of both growing and decaying modes in the gravitational waves. We want the condition 
(assuming AT/Tr is positive; in the case that it is negative the inequalities must be reversed) 



< 



< | AT/Tr | 



(305) 



aF (vr) 

& similarly for ft, 7. In its full form this equation is transcendental even when discussing weak 
waves, but expanding (| 155[) to lowest surviving order in rj, we obtain 



|37C„ M (v% - V 2 e) + 4C a2 ,! ( Vr 3 - Ve 3 ) I < I AT/Tr I 



(306) 



. In a young universe, the times of emission and reception of a light ray are related by t\e 

TjR [Z + l) SO 



37C„i,i(l -(« + !) 1 )»^ + 4C a a, 1 (l-(2 + l) a/: ')»Ja a < I AT/Tr 



\3/2 



(307) 



Let: 



10 -9 be the amplitude of the growing mode C a ^\, so C Q i.i = sgn (C a i t i) 10~ 9 ; 
10~ d be the amplitude of the decaying mode C Q 2,i, so C a %,\ = sgn(C a 2,i) 10 _d ; 
10 _fc be the value of ?7r; 
10" T be the value of \AT/T R \ 



so noting that z ~ 1000 = 10 3 our condition becomes approximately 



S gn(C ahl )10- 2b -^ 2 -sgn(C a2 , 1 )10 3b - d+5 \ < \lQ- T \ 



(308) 



. When the amplitude of the growing mode term dominates, this approximate inequality is 
satisfied by 
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— 2b - g + 3/2 < — T (309) 
; when the decaying mode dominates, the inequality is satisfied by 

36-d + 5<-T (310) 

. WMAP constrains T « 4 (the difference between lowest and highest temperatures is 2ATr/T = 
1.4 x 10~ 4 ) and b > 1. This constrains the growing and decaying modes, when they act on their 
own, to: 



.9>7/2 (311) 
d>12 (312) 

. There exists a third possibility, in which the growing and decaying contributions are, in the 
current epoch, of equal size and opposite sign. For this to be the case, we need 



-2&-.g + 3/2«36-d + 5 (313) 

; this approach relies on the observation of amplification of weak gravitational waves in rapidly- 
changing backgrounds (see section I11.8P . Since b is a free parameter this approximate equation 
can always be satisfied, but we still need to satisfy the constraints of the CMB. 



19 Acceleration in the Bianchi IX universe 

19.1 Order of magnitude estimates for gravitational wave amplitudes 

Meanwhile, consider the tensorial deceleration parameter: 



(an) 



2 



Qo - 2^a - (a + a 2 ) 



■2 



l + 2 ^a + ^a 2 ) (314) 



-l 



& similarly for Q^Q^; this relationship is exact. Evaluating ()225l) gives to lowest surviving order 
in r\ 

AQ (1) . «-70C Q nf7 2 (315) 
^ n, growing al ' 1 1 y ' 

(!) 1^ 
^^n, decaying ~ 2 ^ a2 ir ^ (316) 



With an observed AQu «-lwe can write: 
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10° « sgn (C a i,i) lO 9 / 5 " 26 ^ - sgn (C«a,i) lO 1 ^ 6 (317) 

. In the case of the growing mode dominating we need sgn (C a i,i) = +1 and 9/5+26— g ps 0. This 
forms a system of equations with (|309[) so we have, at the limit of the allowed CMB perturbation, 



J 9/5 — 26- 5 psO , . ,„,„, 

I ' y =^ no solution 318 

1-26-5 + 11/2 « 

; the growing mode cannot, on its own, cause the observed acceleration and be compatible with 
the CMB. For the decaying mode, we need sgn (C a 2,i) = —1 and have 



1 - d + 6 ps I 6 ss -5 J77-IXIO 5 

36 - d + 5 ps -4 ^ I d w -6 ^ 1 C a2 ,i ~ 1 x 10' 



(319) 



which is a nonsense result. Therefore neither the growing or decaying modes, on their own, 
can both cause observed acceleration and preserve the CMB. In the cases of the two modes 
having comparable effect on the metric and opposite sign, though, we can solve (|317p with 
sgn(C Q i i i) = +1, sgn(C a2 ,i) = -1 and 

10° ps - (C qM ) lQ 9 / 5 - 2b -9 + (C a2 ,i) lQ 1 -^ (320) 



] g-dw -56-7/2 , 17 „ „, 52 

{ , =^d~ 777 + 76,gw 26-— (321 

|9/5-26-g«0 10 ,J 10 y ' 

when the growing mode dominates the change in acceleration; this sets estimated limits on the 
parameters (since 6 > 2): 



C a i,i >2 x 10 1 (322) 
C q2 ,i <2 x 10" 16 (323) 



When the decaying mode dominates the change in acceleration 
\g- dps -56-7/2 



=► 5 «-46-5/2,d«6 + l (324) 
1 — a + 6 ps 



which constrains the parameters 



Cod.i >3 x 10 10 (325) 
C a 2,i <1 x 10~ 3 (326) 
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. While the values for the growing mode are far greater than those for what could be called 
"weak" waves (recalling the constraints of section lll.4p . our educated estimate for C a i i in 
the growing-mode dominated regime aligns nicely with the necessary strong-wave growing-mode 
value for AQ\ disregarding the CMB. Therefore we can turn to an analysis in the quasi-isotropic 
regime. 



19.2 Quasi-isotropic, strong growing mode acceleration 

We apply the same reasoning as in the previous section, but we are aware of constraints (from 
[28]) not just on the CMB in the direction of the observed acceleration (which we continue to 
assign as the "a" or e\ direction) but on the CMB in the other two (the "beta" and "gamma" 
directions) : 



AT a /T R + - (ATp/T R + AT y /T R ) =7.1 x 1(T 5 (327) 

2AT (j /T R =1.4 x 1(T 4 (328) 

2AT 1 /T R =1.4 x 10~ 4 (329) 

Q}=- 0.6 (330) 

VR ^3 x 1CT 2 (331) 



In this and all regimes to follow we can also approximate Qp ps Q f = 1/2 to the limit of 
precision given the constraints on ij; Qp will be 1% stronger than Q F A ^ only when rj rj 0.51. 
Between the constraints (|327H330j) and the average over the sky of AT/T R = 0, we have four 
equations with seven unknowns (77, Cq , Cq,Cq, C a 2, C/32, C^). These equations are, explicitly (see 
equations HH [2551 [M|) : 



7.1 x 10" 5 > (ve/Vr) 2 (e o( "* ) - a(,w) - l) 
1.4 x 10~ 4 > (ve/vr) 2 (eP^-fti** - l) 
1.4 x 10~ 4 > (ve/vr) 2 (e 7( " E) - 7(,)R) - l) 

! Qp- tan (gg/2) a' R - tan 2 (gg/2) a" R - tan 2 ( Vr /2) q| 



1 + 2 tan (n R /2) a' R + tan 2 ( m /2) a.% 



(332) 
(333) 
(334) 

(335) 



Trivially, we can see that in the limit of a,/3,~/ —> 00, we must have Q\ Q\ m — > —1: if 
acceleration is driven by growing modes of long-wavelength gravitational waves then in the long 
run the universe asymptotically approaches de Sitter expansion as if driven by a cosmological 
constant, indicating a solution in the regime of quasi- isotropy. 

Consider the quasi-isotropic solution to the growing mode of the Einstein equations, normalized 
as in equations (|215H220j) . In the regime where Cq is sufficiently large that A ^> 1, we can 
approximate 
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« - Jv4 2 (336) 
« <3 ^ 2 (337) 

(an identical argument, with the functions a and (3 transposing their roles, applies for the case 
where Cq < 0). From these terms we can also approximate the next order terms in the series: 



=f»cJ«-^A« (339) 



Approximating equation (|335| to order A 4 ?7 4 we obtain the relationships 



,i {m) ®f - (m-o +m as\ + 

1 - |A 2 ?7| + (is + Tloo) Ai VR 



_9_ 42^2 _l f JT5 9_ 

)_ 

1 + W A2f lR + (lOO ~~ 5e) "^fl 



Q 2 (»W ~ Q-3 lm) - ,,3 „ 9 _ a , / 9 15 \ 44 ^4 




. When Q} (t/r) = -0.6 then Ary fl w 1.5 ± 0.2 (eg* > 1.9), within the limit of applicability of 
the expansion and also in the regime where the infinite series (|212|) converge. Thus, we have 
shown analytically that long-wavelength gravitational waves can explain cosmic acceleration if 
that acceleration is anisotropic. 



We can also make the following qualitative assessments about acceleration. Firstly, its time- 
evolution is non-monotonic. In the a direction, the universe will at first exhibit slightly increased 
deceleration, before starting to accelerate. In the f3 and 7 directions, deceleration will asymptoti- 
cally increase toward infinity but then acceleration will decrease from infinity, quickly converging 
on the strong-field value of Q\ = Q\ = — 1. Acceleration in the a direction begins at Aij « 1.2 and 
the universe accelerates in every direction after Ar\ 1.6; thus the supposition that acceleration 
is a recent phenomenon is supported. 

A universe that is accelerating in every direction is within the region allowed by the model. 
Figure <J2J) illustrates the evolution of the deceleration parameters as a function of time. The 
constraints placed on the decaying mode in (section [T7|) and the upper limit on 77^ show that the 
decaying mode of long- wavelength gravitational waves has not played a significant role in cosmic 
acceleration; in the epoch of last scattering, the deceleration parameter was almost isotropic and 
had a close to Friedmannian value. 



We now turn our attention to the preservation of the CMB. We have three equations in three 
unknowns, taking the lowest term in the decaying mode and the lowest two terms in the growing 
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Figure 2: Deceleration parameter versus time 
Along one direction, the universe at first decelerates, then quickly begins accelerating. Along the 
other two directions, the deceleration parameter goes to infinity before converging from negative 
infinity to the value — 1. The vertical axis of each graph gives Q b a \ the horizontal axis is in units 
of Arj. 
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mode: 



7.1 x l(T a >4 (ve/vr) 



1.4 x io- 4 >4 { VE /mY 



1.4 x 10- 4 >4 (r) E /m) 



c% (v E -V R ) + H^) 2 (v 2 E~r 1 R) 2 + 
+ct (v 4 e ~ Vr) + C a2 ,i (v E 3 ~ Vr 3 ) 

4(VE-VR) + ±(4) 2 ( V l~r lR ) 2 + 

+4(v%-Vr)+C^i(v e 3 -Vr 3 ) 

cUv 2 e~V 2 r) + H^) 2 (v 2 e-Vr) 2 + 
+cj (v E ~ Vr) + C 7 2,i (Ve 3 - Vr 3 ) 



(340) 
(341) 
(342) 



. As 30t]e ~ tjr we can further approximate 



7-1 x 10- >- 



1.4 x 10~ 4 > 
1.4 x 10 -4 > 



1 

225 
1 

225 



20 



20 



A 2 vl 



A 2 vl 



9 | 15 
800 + 224 

9 | 15 
800 + 224 



A 4 ?? 4 + 27000G> 2 , 1 7^ 3 
A\ A R + 27000C 72 ,i7^ 3 



(343) 
(344) 
(345) 



. If we take the inequalities as approximate equivalences and use At\r « 1.5 then this system has 
solutions 



C/32,1 
C Q 2,1 



72,1 



-3 x 10- 7 4 
-9 x 10- 6 4 



which is compatible with the estimates of (section 119. ip . That C a 2,i + C/32,1 + C 7 2,i ^ is 
a consequence of the impossibility of a priori choosing an "unperturbed" temperature against 
which to compare anisotropic CMB fluctuations; the significance of non-linear terms means we 
cannot at the same time have the average over the whole sky of AT/Tr — and have 61 = 0, 
recalling (equation I270p . 

We exhaust almost all the freedom in the system (|332ll335j) in choosing to explain the "Axis 
of Evil" at the same time as acceleration; if this requirement is dropped and we treat CMB 
variations as insignificant then a broad range of solutions opens up. In particular, if the ratio of 
growing mode to decaying mode is approximately equal for all three of a, /?, 7 we always have 
sufficient freedom to choose a 77 that reduces CMB variation to below the level of detectability, at 
the expense of "tuning" the universe to place us as observers in the era when the CMB is nearly 
isotropic. 



Compatibility with an almost-isotropic Hubble flow 



The objection could be raised that the necessity of the universe contracting along two axes 
demands that a large region of the sky be blue-shifted, which would surely have been observed. 
This problem can be made to vanish into statistical noise by the choice of a sufficiently small 77 
as (|223p implies apHu = a' F /ap + a' « 2 (rj^ 1 + c^rj) etc. 
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20 Conclusions 

It is possible for a Bianchi IX universe with initial conditions Cq , Cg , Cq ~ 1 to display the 
acceleration observed in our universe while not only remaining compatible with the observed 
CMB but providing and explanation for potentially meaningful patterns in the CMB, specifically 
the so-called "Axis of Evil" and its associated phenomena such as cold spots. These conditions 
can be attained without additional constraints on the cosmological parameter of Qk , a parameter 
which is subject to further scrutiny and potentially tightening toward the flat universe case of 
VL K = 0. 

The method of combining strong growing modes with linear-order weak decaying modes of cosmo- 
logical gravitational waves is borne out by observational data, which imply a difference of at least 
17 orders of magnitude in amplitude between the growing and decaying modes. In the current 
epoch, decaying modes of cosmological gravitational waves can be neglected entirely. However, 
in the time close to last scattering, these modes may have participated at a strength comparable 
to the growing modes. Furthermore, the action of growing or decaying modes on their own is 
ruled out as an explanation for acceleration as neither on its own can preserve the CMB. 

The question of how the ratio of growing mode to decaying mode is equal along all three principle 
axes of the metric tensor is answered easily if we postulate that cosmological gravitational waves 
present at the earliest moment in time were all in phase (the easiest way to do this is to postulate 
that they consisted of pure growing modes). As the functions a, /3 and 7 would have all crossed 
the boundary from a w = 1/3 medium to a w = medium at the same time, they would thus 
have remained in phase after last scattering, implying equal growing-to-decaying ratios for all 
three functions. As this transition happened in the very young universe (tje ^ 2 x 10 -3 ), the 
decaying mode that exists after last scattering would be very small. 

The nonlinearity of Bianchi IX causes growing modes with initial values of order unity to develop 
exponentially and cause very powerful effects. The structure of the equations also indicates that 
multiple sets of initial conditions can lead to the same set of cosmological parameters. In light of 
the requirement of this model that both strong growing modes and weak but non-zero decaying 
modes of the gravitational waves exist, the possibility that these long- wavelength gravitational 
waves constituted the "pump field" of inflation [103] in the early universe should be explored. 

The model proposed can be tested and falsified by observation of acceleration in areas of the 
sky 90° from the highly-observed field; in areas of the sky away from the currently-observed 
acceleration, we will see either a very large deceleration parameter or a negative one. From the 
analysis of acceleration data in (Part[TI| it is easy to see that, in the current state of observations, 
there are several possible areas of the sky where evidence of a gravitational- wave nature of cosmic 
acceleration could be sitting undetected. 
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21 Directions for future research 

The possibility of explaining cosmic acceleration through a Bianchi IX cosmological model opens 
up numerous possibilities for future research, both theoretical and observational. 

While the difficulties with carrying out a full-sky optical survey of supernovae are understandable, 
experimental verification or falsification of a Bianchi IX model for acceleration requires nearly 
full sky coverage at high z to discover or rule out regions of anisotropy in the acceleration 
field. Infrared astronomy with wide sky coverage, for example WFiRST[65j, presents the best 
possibility for these new observations through traditional astronomy. The Einstein telescope 
provides the tantalizing possibility of independent verification of the properties of dark energy 
through the examination of gravitational radiation. [1 1 1] 

Meanwhile, the available supernova data can be re-examined for signs of acceleration, although 
given the comparatively small datasets in any particular area other than the highly-observed 
field and the equatorial bias in the distribution of the data this re-examination is less likely to 
produce definitive results. Celerier is justified in her criticisms [91] of the assumptions being made 
in proposed models of cosmological acceleration; it is curious that the authors if [89] reasoned, 
with 44 low-z sources, that "poor coverage at low and moderate Galactic latitudes [...] makes it 
practically impossible to distinguish between a peculiar monopole and a quadrupole" but that 
PQ, which shares two authors with [89], does not even mention the possibility of cosmographic 
bias in its smaller sample of high- 2 sources. 

Consideration should be given to the question of why cosmographic bias exists, and whether it 
points to an unexpected privileging of the observer: namely, the fact that modern observatories 
are hosted only in regions of the Earth that can afford to host them. 

Perturbative methods for solving the Einstein equations for weak gravitational waves in Bianchi 
IX can be considered exhausted, having reached the limit of practical utility at quadratic order. 
Further analytic explorations should concentrate on the quasi-isotropic approach. The fact of 
Bianchi IX's easy reduction to a system of non-linear second-order ordinary differential equations 
combined with the divergence of Taylor series describing strong gravitational waves point toward 
either a Fourier-series approach or numerical methods for further analysis; the likelihood of 
chaotic behavior [73] in Bianchi IX, though, merits caution in the selection of initial conditions 
for any simulation. 

Numerical examination of the quasi-isotropic regime should also be pursued for a fuller explo- 
ration of the space allowing for anisotropic acceleration while preserving an almost-isotropic 
cosmic microwave background. The next generation of microwave anisitropy probe should set- 
tle the question of whether the "Axis of Evil" and similar phenomena are genuine artifacts or 
statistical noise; in the meantime, the question of pareidolia in relation to the CMB has not 
been explored and deserves formal examination in order to raise awareness within the scientific 
community of the issue. 

Overall, any theory is only as good as its ability to predict future results. Cosmic acceleration 
needs to be more closely examined, not only for time dependence, but for spatial dependence, 
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before any theory can emerge as preferred. 



22 Implications of the Bianchi IX cosmological model 

Since the discovery of cosmic acceleration, a wide range of scalar theories, ranging from the 
mundane to the exotic, have been put forward to explain the phenomenon. While the fact of 
acceleration, the discovery of which was the logical culmination of the hunt for the "missing mass" 
of the universe above and beyond that provided by dark matter, necessarily implies the slaughter 
of at least one sacred cow, the community of physicists has no consensus over which should be 
sacrificed the most readily. 

Attempts to surrender homogeneity are physically the best-grounded but philosophically the 
most rash. Certainly the idea of a purely homogeneous cosmology is an approximation; but a 
universe which is not on average homogeneous, that is, where the homogeneous regions are rare 
exceptions, is one in which cosmology as a science ceases to be possible. The "Swiss cheese" 
universe has the advantage of making use of a known, exact solution to the Einstein equations 
and at least avoids the exceptionalism of the "Hubble bubble" proposal, but defeats itself on the 
grounds of testability. 

Meanwhile, postulation of exotic states of matter has been done too enthusiastically for the evi- 
dence available. The simple fact of noting that the available data on acceleration was anisotropic 
exposes as irrational exuberance the rush to explain the phenomenon through the medium of a 
substance which has never been seen or even indicated in the laboratory, and whose theoretical 
justification is far beyond testability. The willingness of many to see acceleration as a falsifica- 
tion of the theory of general relativity looks all the more bizarre when counterposed with the 
unwillingness to explore gravitational-wave solutions to the problem. 

The objection could be raised that asserting acceleration to potentially be anisotropic, in the weak 
sense of the word "isotropy", violates the cosmological principle by saying that our telescopes are 
privileged observers, in that our observational field happens to align with an axis of acceleration. 
This is no more so true than the "privilege" hypothesized by, for example, Riess et. al. when 
they assert, from a few dozen data points, that acceleration is a recent phenomenon, and that 
implicitly we are privileged observers in time for taking up cosmology just as the universe has 
begun to exhibit this behavior. While a cosmological constant is the simplest explanation for 
wx = — 1 on mathematical grounds, the lack of physical justification for a non-zero cosmological 
constant puts it in the same class as scalar-field theories. The simple fact is, wx = —1 is, in 
the long run, the natural equation of state for any function which grows faster than the matter- 
driven terms in the background cosmology. The idea of the "Big Rip"[99j, while intellectually 
(and emotionally) intriguing, makes the same mistake in the other direction, privileging observers 
to be alive just as the universe is beginning to tear itself apart. In this sense, a wx = —1 field 
is the best preserver of the cosmological principle, and when the cosmological constant has been 
excluded the simplest explanation for acceleration comes from a tensorial field. 

Similarly, when cosmic flatness is called into question - and it has never been and can never be 
definitively proven, it can only be disproven - the next-simplest model is the closed model. Recall 
that the Bianchi models are distinguished by their symmetries, and of all the Bianchi models with 
Friedmann universes as special cases, Bianchi IX has the most symmetric symmetries, obeying a 
"handedness" rule students learn before their first year of university. The fact of this "handedness" 
- parity - may even provide a neat explanation of the CP violation in particle physics [100) . as 
Grishchuk alluded to[TU]. 
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The least speculative fact revealed by the assessment of available acceleration data is that more 
data is needed, from broader areas of the sky. The anticipated launch of WFIRST is likely to 
prove more momentous for cosmology than the flight of WMAP; WMAP largely reconfirmed 
what we already believed we knew, but WFIRST and SNAP will clearly illustrate how much 
we do not know. We also need techniques to see deeper into the sky and measure the distance- 
redshift relationship further into the past; the standard ladder of baryon acoustic oscillations[101J 
combined with better redshift data from WiggleZ may provide the necessary window. 

That Bianchi IX could in principle contain accelerating regimes was never really in doubt. Nu- 
merical and qualitative analysis has indicated this ever since [10J noted that the vacuum equations 
contained a regular minimum, implying a positive first derivative for the Hubble parameter. The 
character of the acceleration has now been more properly investigated, bringing with it the pos- 
sibility of a purely gravitational explanation for inflation, especially in light of the divergence of 
S constructed only from growing modes in the radiation-dominated universe. An exploration of 
the differences between Bianchi I and Bianchi IX in a universe filled with ultra-relativistic matter 
could make Bianchi IX into a panacea for all the major problems of large-scale cosmology. 

The unwillingness of the perturbed Bianchi IX cosmology to support decaying- mode gravitational 
waves stronger than linear order is puzzling, especially as the BKL universe always has a divergent 
term. The BKL universe, though, never reaches a singularity, and so the divergence of the a 
decaying mode never has time to take effect. Furthermore, the power law contraction along one 
axis could always be explained by a "growing" (non-diverging) function with negative coefficients, 
due to the exponential term in the metric. 

The impact of strong waves on the CMB, meanwhile, also requires deeper explanation. Preser- 
vation of the CMB's apparent anisotropy at first glance appears to require some "tuning", a 
particular growing-decaying ratio which merits deeper questioning; there is also the outstanding 
matter of why we happen to live in one of the few periods of time when the CMB appears nearly 
isotropic. Clever examination of the symmetries of Bianchi IX may reveal a more satisfying an- 
swer, although the ability of Bianchi IX to explain CMB anomalies is one of its most satisfying 
features. 

Most fundamentally, the biggest impact of the Bianchi IX theory of cosmic acceleration is the 
expansion of the cosmologist's parameter space. While in scalar models the only parameter truly 
open for discussion is the function describing the equation of state of dark energy, the gravi- 
tational waves of the Bianchi IX universe have four degrees of freedom; while the strength a 
non-zero cosmological constant has some theoretical justification in fundamental physics inde- 
pendent of large-scale cosmology, there is no immediately apparent reason why the gravitational 
waves in Bianchi IX should have any particular amplitude. As always in cosmology, we need 
more information than we have. 
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Table 3: Supernova observations used in analysis of acceleration 





Right ascension, J2000 


Declination, J2000 


Riess 1998 supernovae:T9(N 


SN1994U 


13:04:56 


-6:3:39 


SN1997bp 


12:46:54 


-10:21:27 


SN1996V 


11:21:31 


2:48:40 


SN1994C 


07:56:40 


44° 52' 19" 


SN1995M 


09:38:42 


-11:39:52 


SN1995ae 


23:16:56 


-1:55:24 


SN1994B 


08:20:41 


15:43:49 


SN1995ao 


02:57:31 


-0:18:40 


SN1995ap 


03:12:28 


0:41:43 


SN1996R 


11:16:10 


0:11:39 


SN1996T 


10:05:28 


-6:32:36 


SN1997I 


04:59:37 


-2:50:58 


SN1997ap 


13:47:10 


2:23:57 


SDSS-II SNIa observations:^ 


(Corner 1) 


20:00:00 


1:15:00 


(Corner 2) 


20:00:00 


-1:15:00 


(Corner 3) 


04:00:00 


1:15:00 


(Corner 4) 


04:00:00 


-1:15:00 


ESSENCE windows:^ 


waal 


23:29:52.92 


-08:38:59.7 


waa2 


23:27:27.02 


-08:38:59.7 


waa3 


23:25:01.12 


-08:38:59.7 


waa5 


23:27:27.02 


-09:14:59.7 


waa6 


23:25:01.12 


-09:14:59.7 


waa7 


23:30:01.20 


-09:44:55.9 


waa8 


23:27:27.02 


-09:50:59.7 


waa9 


23:25:01.12 


-09:50:59.7 


wbbl 


01:14:24.46 


00:51:42.9 


wbb3 


01:09:36.40 


00:46:43.3 


wbb4 


01:14:24.46 


00:15:42.9 


wbb5 


01:12:00.46 


00:15:42.9 


wbb6 


01:09:00.16 


00:10:43.3 


wbb7 


01:14:24.46 


-00:20:17.1 


wbb8 


01:12:00.46 


-00:20:17.1 


wbb9 


01:09:36.40 


-00:25:16.7 


wccl 


02:10:00.90 


-03:45:00.0 


wcc2 


02:07:40.60 


-03:45:00.0 


wcc3 


02:05:20.30 


-03:45:00.0 


wcc4 


02:10:01.20 


-04:20:00.0 


wcc5 


02:07:40.80 


-04:20:00.0 


wcc7 


02:10:01.55 


-04:55:00.0 


wcc8 


02:07:41.03 


-04:55:00.0 


wcc9 


02:05:20.52 


-04:55:00.0 


wdd2 


02:31:00.25 


-07:48:17.3 


wdd3 


02:28:36.25 


-07:48:17.3 
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wdd4 


02:34:30.35 


-08 


19:18.2 


wdd5 


02:31:00.25 


-08 


24:17.3 


wdd6 


02:28:36.25 


-08 


24:17.3 


wdd7 


02:33:24.25 


-08 


55:18.2 


wdd8 


02:31:00.25 


-09 


00:17.3 


wdd9 


02:28:36.25 


-09 


00:17.3 


HST supernovae:\H\ 


Dbr UoiJU 


no. 01 • A O Cifld 
Uz.Zl.4z.U00 


HQ -01 • 10 

— Uo.zl.oo.lz 


Dbr UOrio 


1 A .QH 1 Af\ 
14.o4.oU.14U 


o A -Oft- K*7 QH 

o4:zd.o < .oil 


OV^r UOrvU 


14.oo.Uo.oD0 


Q/1 ■ 1 A -1 HQ 
o4.14.1o.UO 


SCP06K18 


14:38:10.665 


34:12:47.19 


SCP06R12 


02:23:00.083 


-04:36:03.05 


SCP06U4 


23:45:29.430 


-36:32:45.75 


SCP06C1 


12:29:33.013 


01:51:36.67 


SCP06F12 


14:32:28.749 


33:32:10.05 


SCP05D6 


02:21:46.484 


-03:22:56.18 


SCP06G4 


14:29:18.744 


34:38:37.39 


SCP06A4 


22:16:01.078 


-17:37:22.10 


SCP06C0 


12:29:25.655 


01:50:56.59 


SCP06G3 


14:29:28.430 


34:37:23.15 


SCP06H3 


14:34:28.879 


34:27:26.62 


SCP06N33 


02:20:57.699 


-03:33:23.98 


SCP05P1 


03:37:50.352 


-28:43:02.67 


SCP05P9 


03:37:44.513 


-28:43:54.58 


SCP06X26 


09:10:37.888 


54:22:29.06 


SCP06Z5 


22:35:24.967 


-25:57:09.61 


Riess "gold" dataset:j56^ [54V 


Window 1 


03:32:30 


-27:46:50:00 


Window 2 


12:37:00 


62:10:00 
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